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ABSTRACT 


This document presents the results of an extensive survey 
and comparative evaluation of current atmosphere and wind 
models for inclusion in the Langley Atmospheric Information 
Retrieval System (LAIRS) . It includes recommended models 
for use in LAIRS, estimated accuracies for the recommended 
models, and functional specifications for the development 
of LAIRS. 
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SECTION 1 - INTRODUCTION 


Postflight determination of aerodynamic force coefficients 
for the Shuttle Orbiter via the Aerodynamic Coefficient 
Measurement Experiment (ACME) will require independent 
knowledge of the characteristics of the atmosphere along the 
Shuttle reentry path. 

As a part of the planning for the ACME Program, the NASA 
Langley Research Center is making arrangements with the 
Shuttle Project to obtain certain atmosphere measurements 
near the time and path of each reentry flight. These 
measurement data, taken from balloons, sounding rockets, 
aircraft, and ground stations, will be processed by LaRC 
to determine the needed atmospheric velocity, pressure, 
temperature, and density along the Shuttle trajectory-. 
Knowledge of the atmosphere's local velocity vector not 
only will permit the determination of the relative speed, 

V, but, together with the vehicle attitude, will also make 
possible the computation of the aerodynamic angle of 
attack, CL, and the sideslip angle, /3 . Knowledge of 
thermodynamic properties other than density will define 
parameters such as the local speed of sound and the viscos- 
ity. These parameters will in turn determine the flight 
Mach number and Reynolds number for correlation of the 
aerodynamic coefficients with other experimental data. 

In support of this activity, an atmospheric data processor, 
the Langley Atmospheric Information Retrieval System (LAIRS) , 
will be developed for use in the ACME Project. In prepara- 
tion for the development of LAIRS, Computer Sciences 
Corporation (CSC) has conducted a survey and comparative 
evaluation of current lower atmosphere, upper atmosphere, 
and wind models; recommended models for use in LAIRS; 
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estimated the expected accuracies of the recommended models, 
and generated functional specifications for LAIRS. This 
report contains the results of these studies . 

Section 2 discusses the various atmosphere and wind models 
considered and the models recommended for use in LAIRS. 
Section 3 discusses diurnal, semidiurnal, and latitudinal 
variations and presents procedures for fitting temperature, 
density, pressure, and wind profiles. It also includes 
a discussion of estimated model accuracies. Section 4 con- 
sists of the functional specifications for LAIRS. 

The authors wish to acknowledge the direction and assistance 
of Gloria P. White of CSC, who contributed many helpful sug- 
gestions and comments, and the overall guidance of the work 
by Joseph M. Price of the Langley Research Center. 
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SECTION 2 


EVALUATION OF ATMOSPHERIC MODELS 


Currently available atmospheric models are reviewed in this 
section. Atmospheric models that are appropriate for the 
lower part of the atmosphere are discussed in Section 2.1 
and those for the upper part of the atmosphere are discussed 
in Section 2.2, while Section 2.3 deals with wind modeling. 

In each case, recommendations are made for specific models 
to be implemented in the Langley Atmospheric Information 
Retrieval System (LAIRS) . 

2.1 LOWER ATMOSPHERE MODELS 

The Earth's atmosphere from ground level up to approximately 
100 kilometers is referred to as the lower atmosphere in this 
report. This portion of the atmosphere consists of several 
layers with different physical properties. The troposphere, 
which is the region nearest the Earth's surface, exhibits a 
more or less uniform decrease of temperature with altitude 
and is bounded by the next layer, the stratosphere, at an 
altitude of about 10 kilometers. The troposphere is the 
domain of weather and is in convective equilibrium with the 
Sun-warmed surface of the Earth. The stratosphere has about 
the same vertical dimension as the troposphere and has a 
nominally constant temperature. The stratosphere is known 
to have mostly horizontal winds. The mesosphere lies above 
the stratosphere. It is the region containing the first 
temperature maximum and is bounded from above by the major 
temperature minimum, which is found near an altitude of 
85 kilometers. The thermosphere covers the entire atmo- 
sphere above the mesosphere and is characterized as a region 
of rising temperature. A typical temperature-altitude pro- 
file is shown in Figure 2-1. 

If the atmosphere were macroscopically motionless for a 
sufficiently long time, it would show a decrease in mean 
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molecular weight with altitude, with the lighter species 
having separated from the heavier species by the process of 
molecular diffusion. Such a process is not observed in the 
Earth's atmosphere until an altitude of about 100 kilometers 
is reached; up to this altitude, the mean molecular weight 
remains approximately constant. This comes about because 
the Earth's atmosphere is not motionless on a macroscopic 
scale, but is continually mixed by atmospheric turbulence 
and large-scale circulation. The beginning of diffusive 
separation at altitudes near 100 kilometers does not mean 
that there are no winds to cause turbulence and mixing above 
these altitudes, but rather that at such low densities molec 
ular diffusion rates are much greater than the wind-induced 
'mixing rates. 

The lower atmosphere is thus assumed to have a uniform com- 
position, in the sense of having a constant mean molecular 
weight. It is also assumed that the air behaves as if it 
were a perfect gas obeying the perfect gas law: 


p = 

where 


P 

T 

M 


total pressure 
total density 
absolute temperature 
mean molecular weight 


R* = universal gas constant 

Within the height region of complete mixing, the atmosphere 
is further assumed to satisfy the hydrostatic equation, i.e. 


cLp =-pgdz 
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where 


acceleration of gravity 
geometric altitude 


g 

Z 


Equations (2-1) and (2-2) can be combined to obtain the 
integral form of the barometric equation: 

P ^dZ 

2 . 

where 

p = pressure at Z 

Pq = pressure at Z^ 

H = R*T/Mg 

P 

The quantity is referred to as the pressure scale height. 
The total atmospheric density j> can be obtained from Equa- 
tions (2-1) and (2-3): 



where 

f. = M 

The density scale height, Ej> , can similarly be defined through 
the following equation: 






(2-5) 
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The density scale height and the pressure scale height are 
related through the following equation: 




:e 1 il 

Mg dz 




( 2 - 6 ) 


where M is treated as a constant. 

The following atmosphere models that are currently available 
were developed using the assumptions discussed above: 

1. ARDC Model Atmosphere 1959 (ARDC 1959) 

2. U.S. Standard Atmosphere 1962 (USSA 1962) 

3. Cospar International Reference Atmosphere 1965 
(CIRA 1965) 

4. Cospar International Reference Atmosphere 1972 
(CIRA 1972) 

5. Goddard Institute for Space Studies 1974 Atmosphere 
(GISS 1974) 

6 . The Edwards Air Force Base Reference Atmosphere 
(ERA 1975) 

7. U.S. Standard Atmosphere 1976 (USSA 1976) 

These models are discussed in the following section. 

2.1.1 EVALUATION OF AVAILABLE MODELS 

Atmospheric models are often described in terms of the 
molecular-scale temperature, T^^, and the geopotential height, 
H. The relation between Tj^ and T, the ordinary kinetic 
temperature, is given by 
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where 


M = M(Z) = mean molecular weight of the atmosphere 
at height Z . 


Mq = M(0) 

Since M remains practically unchanged up to about 100 kilo- 
meters, and T are essentially the same for the lower 
atmosphere and begin to deviate from each other above 
100 kilometers. The geopotential height, H, is defined in 
termS; of the geometric height, Z, as follows: 


„ . Cm ,, 


( 2 - 8 : 


where 


g(s) 


H 


go = 9 ( 0 ) 


geometric height measured from mean 
sea level 

acceleration of gravity at Z 
geopotential height 
acceleration of gravity at Z=0 


The inverse-square law of gravitation provides an expression 
for g(Z) with sufficient accuracy for most model-atmosphere 
computations as follows : 


3(^)=3o 




Toi-Z 



(2-9) 
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where Tq is the effective radius of the Earth at a particular 
latitude. Using Equation (2-9) in Equation (2-8) , an approx- 
imate relation between H and Z can be obtained: 


H 


n + Z. 


( 2 - 10 ) 


The major advantage of using instead of T is that it com- 
bines the variable portion of M with the variable T into a 
single new variable. Similarly, combining the variable por- 
tion of g with Z forms the new variable H. 

Among the available lower atmosphere models, the ARDC 1959 
(Reference 1) , USSA 1962 (Reference 2) , and USSA 1976 (Ref- 
erence 3) models were formulated in terms of the molecular- 
scale temperature and the geopotential altitude . 

2. 1.1.1 The ARDC 1959 Model 

In the ARDC 1959 model, the variation of the molecular-scale 
temperature is defined as a series of connected segments that 
are linear in geopotential height from the ground up to 700 
kilometers. The general form of each linear segment is: 



+ (-M ( ''I 


(2-U) 


where T., , 
M,b 


L' 


M 


molecular- scale temperature at the 
geopotential altitude 

temperature gradient with respect to H 
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2. 1.1. 2 The U.S. Standard Atmosphere 1962 

The USSA 1962 model/ on the other hand, uses Equation (2-11) 

only up to a height of 90 kilometers. Above 90 kilometers, 

T is expressed in geometric height Z rather than H. Thus, 

M 

for Z greater than 90 kilometers, 

=T^,b 

where L denotes the temperature gradient with respect to 
M 

geometric height Z . 

2. 1.1. 3 The U.S. Standard Atmosphere 1976 

The structure of the USSA 1976 model is somewhat more involved. 
The temperature-altitude profile of USSA 1976 was defined as 
follows : 

1. From Z = 0toZ = 86 kilometers, Tj^ is expressed as 

a series of connected segments 
linear in H. 

2. From Z = 86 kilometers to Z = 1000 kilometers, the 

kinetic temperature T and the 
geometric height Z are used to 
define the temperature-altitude 
profile. 

a. From Z = 86 kilometers to Z = 91 kilometers, 

T is constant at 186.8673° K. 

b. From Z = 91 kilometers to Z = 110 kilometers, 

T(Z) is fitted to an ellipse. 

c. From Z = 110 kilometers to Z = 120 kilometers, 

T(Z) is linear in Z with a posi- 
tive gradient. 
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d. From Z = 120 kilometers to Z = 1000 kilometers, 

T(Z) is an exponential function. 

These first three models (ARDC 1959, USSA 1962, and USSA 1976) 
represent the year-round average atmosphere in the mid- latitude 
region over North America. 

2. 1.1. 4 Cospar International Reference Atmosphere 1965 
In CIRA 1965 (Reference 4), an attempt was made to develop: 

1. A set of atmospheres for the region from 30 kilo- 
meters to 80 kilometers for various latitudes and 
times of the year 

2. A set of atmospheres for the region from 120 kilo- 
meters to 800 kilometers for various solar fluxes 
and times of day 

3 . A Mean Reference Atmosphere for the altitude range 
from 30 kilometers to 300 kilometers 

The Mean Reference Atmosphere was defined as follows : 

1. Between 30 kilometers and 80 kilometers, the model 
represents the mean conditions throughout the year 
for latitudes near 30 degrees north. 

2. Between 120 kilometers and 300 kilometers, the 
model represents mean conditions over a 24-hour 
period for a 10.7-centimeter wavelength solar flux, 

Fio ~j t of 150 X 10 w/m^/Hz. 

3. Between 80 kilometers and 120 kilometers, a model 
was developed that provides a smooth connection 
between the two parts of the mean atmosphere . 

Data published by the U. S. Meteorological Rocket Network 
and further data from grenade firings were used in preparing 
the mean and monthly reference atmospheres for CIRA 1965. 

Most of these data were collected during daytime and at 
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latitudes near 30 degrees north. This model gives not only 
temperatures, pressures, and densities, but also winds. 
Eurthermore, the latitudinal variations of various atmospheric 
properties are also included in the model. The tabulations 
of the models are at intervals of 5 kilometers in height, 

10 degrees in latitude, and 1 month in time. In determining 
the variation of temperature with latitude, use was made of 
the wind model in conjunction with the theinnal wind equation, 
especially when insufficient temperature data was available. 

2. 1.1. 5 Cospar International Reference Atmosphere 1972 

The CIRA 1972 model (Reference 5) follows essentially the 
same formulation as CIRA 1965. A few notable differences are 
the following; 

1. The lower part of the atmosphere covers the region 
from 25 kilometers to 110 kilometers. 

2. With more data available, the temperature model of 
CIRA 1972 was obtained without using the thermal 
wind equation. 

3. Additional data were collected and analyzed for 
the diurnal and semidiurnal variations associated 
with atmospheric phenomena. 

The model temperatures of CIRA 1972 were determined on the 
basis of observed temperatures from the relationship 

"7T (9) ^■'T: (6) -t a -t b 9 -r c 

new 


where Q is the latitude in the range 0 degrees to 90 
degrees and a, b, and c were obtained using the method of 
least squares. This procedure was repeated with T^^^ as the 
next T^^ until convergence was obtained. Five-point 
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smoothing of in height, latitude, and season was then 

carried out until a smooth model was obtained. 

The pressure p and the density ^ were calculated from the 
temperature T and the 30-kilometer pressure p^ : 

P = (2-13 


f = M/t*T 


(2-14) 


2. 1.1. 6 The GISS 1974 Model 

The GISS 1974 model (Reference 6) is a global atmospheric 
circulation model developed at the Goddard Institute for 
Space Studies (GISS) . This model is a nine-level primitive- 
equation model; the fundamental equations are the equation 
of motion, the equation of continuity, the equation of state, 
the first law of thermodynamics, the hydrostatic equation, 
and a conservation equation for water vapor. These equations 
are numerically integrated with appropriate sets of initial 
and boundary conditions. The lower boundary is the Earth's 
surface and the upper boundary lies at an altitude of about 
30 kilometers, which is determined by the upper boundary 
pressure of 10 millibars. At the lower boundary, the model 
includes a realistic distribution of continents, oceans, 
and topography. The model also calculates cloud and water 
vapor fields and uses them in the calculations of energy 
transfer by solar and terrestial radiation. 
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2. 1.1. 7 ERA Model 


The ERA 1975 model (Reference 7) is a mean reference atmos- 
phere based on temperature observations made at Edwards Air 
Force Base, California, between the years 1953 and 1967 and 
on previously modeled temperatures, pressures, and densities 
for Vandenberg Air Force Base, California. The model ex- 
tends from the surface to an altitude of 700 kilometers. 
Above 81.75 kilometers, it is identical to the USSA 1962 
atmosphere model. Below 81.75 kilometers, the temperature 
profile was divided into six altitude segments, and a fifth- 
order polynomial was fitted to each segment. Pressures and 
densities were derived from the temperature profile below 
3.25 kilometers and were taken directly from the Vandenberg 
model between 3.25 and 81.75 kilometers. These pressure 
and density profiles were also broken into six segments, and 
fifth-order polynomials were fitted to the logarithms of the 
quantities in each segment. The values of the polynomial 
coefficients, along with a set of profiles calculated from 
them, is given in Reference 7. 


2-12 



2.1.2 RECOMMENDED MODEL FOR THE LOWER ATMOSPHERE 


The Aerodynamic Coefficient Measurement Experiment (ACME) 
(Reference 8) requires the instantaneous values of atmos- 
pheric mass density and wind velocities throughout the Shuttle 
altitude range during reentry. Among the lower atmospheric 
models studied, the CIRA 1972 model best meets these require- 
ments. The lower atmospheric part of the CIRA 1972 model 
gives tabulated atmospheric parameters, i.e., temperature, 
pressure, density, and wind velocities at intervals of 5 kilo- 
meters in height, 10 degrees in latitude, and 1 month in time. 
Furthermore, the meteorological data used in CIRA 1972 has a 
high concentration near a latitude of 30 degrees north, which 
is in the neighborhood of the Shuttle reentry path. 

The lowermost 25 kilometers of the atmosphere are not covered 
by the CIRA 1972 model. For this region, the ERA 1975 model, 
which was designed with the Shuttle reentry in mind, is recom- 
mended. Though no monthly variation is included in this 
model, a summary of the data^ from which the model was con- 
structed is available. These data, in monthly form, should 
allow the amplitude and phase of the yearly variation to be 
easily determined. 

The combination of CIRA 1972 and ERA 1975 will provide a 
general default lower atmospheric model for latitude 30 de- 
grees north during the major portion of the Shuttle reentry 
and a more specific default model during the final descent 
phase into Edwards Air Force Base. The two models merge 
well at 25 kilometers; for example, the difference between 
the ERA temperature and the average of the twelve CIRA tem- 
peratures at that height is only 0.21° K. 

The CIRA 1972 model can easily be replaced by a newer CIRA 
model when it becomes available (Reference 9) . 

^Private communication from J. M. Price of the NASA Langley 
Research Center. 
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The diurnal and semidiurnal effects can be incorporated in 
the following approximate manner. Let f(Z,t) be an atmos- 
pheric parameter that is a function of altitude and local 
solar time, where f(Z,t) is given approximately by (Ref- 
erences 5 and 10) : 

f('Z.t) = l(Z)t/l,(Z)a>s^[t-t,(z)J , 2 - 15 ) 

t /li(z)cosi [t-tzCZ)] 

where 

Aq(Z) = diurnal average of f(Z,t) 

A^(Z) = diurnal amplitude at Z 

A^(Z) = semidiurnal amplitude at Z 

t^(Z) = time of the maximum diurnal effect at Z 

t 2 (Z) = time of the maximiam semidiurnal effect at Z 

The function f(Z,t) can be written in the following form: 

= ^^(Z)-Scz.t) (2-16) 

where 

S(Z,t) = 1 t1?,(z)Cos^[t-t,i‘z)]tT?^(z)cosf[f.tjz)J 

Ki(2)= 
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Aq(Z) can be considered as altitude profiles of respective 
atmospheric parameters tabulated in CIRA 1972. However, at 
altitudes below 50 kilometers, the CIRA temperature, pres- 
sure, and density are biased towards local noon, and a con- 
venient form of f(z,t) for this case is 


f ( z , t ) A f 'S( ) 


(2-18) 


where 

SJ‘Z.)=S('Z.,IZ) 

(2-19) 

= 1- R,(Zjcos(it,«)) 


Thus, below 50 kilometers, tabulated CIRA 1972 atmospheric 
parameters can be used in place of A^ 2 ^^) Equation (2-18) . 

In the region below 25 kilometers, diurnal and semidiurnal 
variations in all parameters are quite small, being generally 
less than 1 percent of the average value of the parameters 
except in the lowest kilometer or two of the atmosphere. 

Thus, it will probably be unnecessary to include diurnal 
and semidiurnal effects in the bottom segment of the model 
atmosphere. It will also be relatively easy to include a 
limited degree of latitude dependence in the lower atmos- 
pheric model. This can be accomplished by tabulating sets 
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of latitude gradients for the. various parameters at certain 
altitudes for various times of the year. These gradients 
can be determined from the 1972 CIRA profiles for altitudes 
above 25 kilometers and from the ARDC Handbook of Geophysics 
(Reference 11) for altitudes below 25 kilometers. Clearly, 
this method allows for accurate calculation of latitude 
dependence only within a relatively narrow region (e.g., 

+10 degrees) centered on 30 degrees north latitude. 

2.2 UPPER ATMOSPHERE MODELS 

The upper atmosphere starts at altitudes near 100 kilometers 
and covers the entire thermosphere, which is characterized 
as a region of rising temperature and diffusive equilibrium. 
As shown in Figure 2-1, the neutral gas temperature has a 
profile that increases from a minimum of about 180°K near 
an altitude of 90 kilometers to a constant value above 300 
kilometers that is called the exospheric temperature. Above 
300 kilometers, the air density is low and molecular thermal 
conduction is rapid enough to eliminate any vertical temper- 
ature gradient, keeping the neutral gas temperature constant 
with altitude. The global mean neutral gas exospheric temp- 
erature increases with the solar extreme ultraviolet radia- 
tion and ranges from several hundred degrees to a few thou- 
sand degrees. The atmosphere below 90 kilometers is almost 
entirely in the molecular state and consists of approximately 
78 percent molecular nitrogen, 21 percent molecular oxygen, 
and 1 percent other minor gas constituents (Reference 3) . 
Above 90 kilometers, in the thermosphere, intense solar 
ultraviolet radiation can break the molecular bonds, creating 
atomic species. The photo-dissociation of molecular oxygen 
into atomic oxygen is particularly effective in the thermos- 
phere. This changes the compositional structure of the 
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atmosphere from a molecular mixture in the lower thermosphere 
to a mixture rich in atomic oxygen in the upper thermosphere. 

There are several types of variation recognized in the 
thermosphere. They can be classified as follows (References 
12 and 13) ; 

1. Variations with solar activity 

2. Diurnal variation 

3. Variations with geomagnetic activity 

4 . Semiannual variation 

5. Seasonal-latitudinal variations 

6. Rapid density fluctuations, probably connected with 
gravity waves . 

Currently available upper atmosphere models have attempted 
to take into consideration as many of these variations as 
possible. Among these available models, the following 
models were chosen for evaluation: 

• Harris-Priester 1962 model 

• Jacchia-Roberts 1971 model 

• Jacchia 1977 model 

• MSIS 1977 model 

• Russian models 

In these models, variations of atmospheric parameters with 
solar activity were formulated in terms of the solar 10.7- 
centimeter radio flux, which is known to be closely related 
to the solar extreme ultraviolet (EUV) radiation, one of the 
main heat sources for the upper atmosphere. The solar 10.7- 
centimeter radio flux can be measured from the ground, while 
the solar EUV radiation cannot be measured at ground level 
because of strong absorption in the upper atmosphere. This 
flux, denoted as FlO.7, gives a measure of the 10.7-centimeter 
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. -22 2 

flux xn units of 10 w/m /Hz. F^q ^ values range from 
about 70 at minimum solar activity to about 250 at maximum 
solar activity. 

2.2.1 EVALUATION OF AVAILABLE MODELS 


2. 2. 1.1 Harris-Priester 1962 Model 


In 1962, Harris and Priester (Reference 10) developed a 
time-dependent model of the upper atmosphere using the 
energy equation and the barometric equation. This model is 
one-dimensional in the sense that the only space coordinate 
considered is the height Z . The second independent param- 
eter in the model is the local solar time, which is equiva- 
lent to the longitude. Latitudinal dependence is completely 
disregarded. The procedure is based upon a perturbation 
approach suggested by Eckart (Reference 14) in which the 
temperature profile at the time t+At is obtained from the 
energy-temperature equation using the known temperature pro- 
file, atmospheric density, and energy input at the time t. 
Harris and Priester used two energy sources, one due to the 
extreme ultraviolet radiation of the Sun and the other to 
account for any other possible sources of energy. 


The second energy source was assumed to have approximately 
the same daily integrated intensity and the same height 
dependence as the EUV source . The time of the diurnal maxi- 
mum of the second source was chosen such that the calculated 
temperature profiles and density variations were in agreement 
with those obtained from satellite drag measurements. There 
are two variations of this original Harris-Priester approach. 
One is the Modified Harris-Priester Model used in the Goddard 
Trajectory Determination System (GTDS) (Reference 15) and 
the other is the Asymmetric Modified Harris-Priester Model 
introduced by Dowd and Tapley (Reference 16) . The Modified 
Harris-Priester Model uses tabulated daytime maximum and 
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nighttime minimum density values at various altitudes and 

^10 7 These tables were constructed from a similar, 

but more extensive, tabulation given in CIRA. 1965 ’(Reference 4) . 

Using the maximum and minimum densities, p (h) and p . (h) , 

^ ' ^max '^min ' 

at an altitude h, the density at h is given by 


( 2 ^ oa ^ 8 ) 

where m is the angle between the position vector of the 
point where the density is desired and the apex of the 
diurnal bulge. This procedure yields a density distribution 
that is symmetric with respect to the apex of the diurnal 
bulge. However, it is known (Reference 17) that the observed 
diurnal variation is not symmetric. To incorporate this 
asymmetric diurnal variation in the Harris-Priester model, 

Dowd and Tapley (Reference 16) proposed the Asymmetric 
Modified Harris-Priester Model. In this model, the asymmetric 
diurnal formulation developed for temperature profiles by 
Jacchia (Reference 13) was used to compute a diurnally asym- 
metric density distribution. The density at h is given by 

= Tn ^ CCS” (4 ) 
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where 


= L i-Qcos">/] 

^ ~ i 1 0 + 5 1 

7 =ile-Jl 

T = H-tyS+X(H + 7} f-7c <T < 7t) 

In these equations, H is the local solar time, Q is the 
latitude of the subsatellite point, and ^ is the solar 
declination. The values of the parameters appearing in 
Equations (2-22), taken from Jacchia's temperature equation 
(Reference 13) , are; 


m = 2 . 2 
3 = -37° 

X = 6° 

Y = 43° 

The two versions of the modified Harris-Priester model 
described above explicitly include in the computation of the 
atmospheric density only the first two variations listed in 
Section 2.2, namely, the effects of the solar 10.7-centimeter 
flux level and the diurnal variation. 

2. 2. 1.2 Jacchia-Roberts Model 

The Jacchia 1971 Model (Reference 12) describes an upper 
atmosphere above 90 kilometers. This model is based upon 
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a family of vertical temperature profiles that were con- 
structed from observed data. The computation of the atmos- 
pheric density involves solving a pair of differential 
equations using an appropriate temperature profile. The 
two. equations are the diffusion equation for the number 
density of individual species above an altitude of 100 kilom- 
eters and the barometric equation for the mean density of 
the entire constituents of the atmosphere between 90 kilom- 
eters and 100 kilometers. The Jacchia-Roberts model imple- 
mented in GTDS (Reference 15) uses analytical solutions to 
these differential equations given by Roberts (Reference 18) . 
The Jacchia-Roberts model includes the effects associated 
with the solar 10.7-centimeter flux, the diurnal variation, 
and geomagnetic activity in the computation of temperature 
profiles. The semiannual and seasonal-latitudinal variations 
were incorporated into the final density calculation. 

Each temperature profile is labelled by the so-called exos- 
pheric temperature, T^ , which sets the upper boundary value 
for the temperature. T^ is computed in three steps in the 
Jacchia-Roberts model. First, from the solar 10.7-centimeter 
flux level, the nighttime minimum exospheric temperature, T^, 
is obtained 

'T, = 37 ?” + ir 3 , 2 - 23 , 

The second step involves computing the exospheric 

temperature which includes the diurnal variation, from T^^ : 

T = ~T t f T -T ) (2-24) 

J>h/L U ^ ^ '2) '/V j _ 
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where 


'l^=7^(li-0.3 Cos"'rj ) 

(2-25) 

Tc ( 1 •+ 0.3 ) 

Various parameters appearing in these equations have already- 
been defined in Equation (2-22) . Finally, in the third step, 
the contribution to from the geomagnetic activity effect 

is computed , i . e . , 

=Z8X (Z>200Km) 

. ( 2 - 26 ) 

= Kf to:ozeXz<zooKm) 


where is the 3-hour geomagnetic planetary index. An aver- 
age time lag of 6.7 hours was used for the time lag associated 
with the atmospheric variation resulting from the geomagnetic 
disturbance . 

The total exospheric temperature is given by 


The value of T^^j uniquely determines the temperature profile 
to be used in the computation of the density, which will be 
further corrected for the semiannual and seasonal-latitudinal 
variations. See Reference 13 for details of these corrections. 
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2. 2. 1.3 The Jacchia 1977 Model 


The general structure of Jacchia 's 1977 model (Reference 19) 
is identical to that of his 1971 model. Temperature profiles 
obtained using the 10.7-centimeter solar flux, the diurnal 
variation, and the geomagnetic activity effect are used to 
compute the atmospheric density from the diffusive equations 
above 100 kilometers and the barometric equation is used 
below 100 kilometers. The differences between the two models 
are primarily in the details of the calculation of the temper- 
ature profiles. 

In Jacchia 1977, the arithmetic mean of the nighttime mini- 
mum and the daytime maximum exospheric temperature, T;^, is 
related to the 10.7-centimeter solar flux index ^ as 
follows : 

\ = 5 . -f 101. 8 1 =^^ (2-28) 

where ^ is a weighted mean of F^^g the weight being a 

Gaussian function of time centered on the current time. 

The diurnal variation of temperature was formulated for each 
constituent of the atmosphere in terms of T^, the declination 
of the Sun, and the hour angle of the Sun. The diurnal temp- 
erature curve has a somewhat more complicated shape than that 
of the 1971 model. The total density maximum and minimum 
occur at 16.8 hours and 5.4 hours local solar time, respec- 
tively. These values are to be compared with 14 hours and 
3 hours local solar time in the 1971 model. 

Also, a much more sophisticated approach is used in computing 
the geomagnetic activity effect in the 1977 model. The 
latitude-dependent time lag is empirically derived with the 
minimum time lag of 2.4 hours at the geomagnetic poles and 
the maximum time lag of 7.2 hours at the equator. Further- 
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more, the exospheric temperature increase due to the geomag- 
netic activity effect is assumed to accompany a change in the 
height of the homopause, which is the point where the atmos- 
phere ceases to be a homogeneous mixture of its various com- 
ponents. The change in the height of the homosphere brings 
about a change in density due to the equatorial wave, which 
is the density pileup in the equatorial regions as a conse- 
quence of convection toward the equator. The semiannual and 
seasonal-latitudinal variations have been taken into account 
in a way similar to that in the 1971 model. 

2. 2. 1.4 The MSIS 19 7 7 Model 

The MSIS model (Reference 20) is a purely empirical model 
based upon a large amount of mass spectrometer and incoherent 
scatter measurement data. These data were provided by mass 
spectrometers on five satellites (Orbiting Geophysical 
Observatory-6 (OGO-6) , San Marco-3, Aeronomy Satellite-1 
(AEROS-1) , Atmosphere Explorer-2 (AE-2) , and AE-3) and four 
incoherent-scatter ground stations (Arecibo, Jicamarca, 
Millstone, and St. Santin) . The exospheric temperature and 
other quantities were represented by expansions in terms of 
spherical harmonics. The lower boundary of the MSIS model 
is at an altitude of 120 kilometers. 

The model coefficients are determined in a series of steps 
that consider independently (1) 120-kilometer temperatures, 

(2) 120 -kilometer densities, (3) exospheric temperatures and 
temperature gradient parameters under quiet magnetic condi- 
tions, and (4) high magnetic activity. 

The MSIS model includes the 10.7-centimeter solar flux effect, 
the geomagnetic activity effect, the diurnal and semidiurnal 
variations, the annual and semiannual variations, and 
seasonal-latitudinal variations. 
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2. 2. 1.5 The Russian Model 


The Russian atmospheric model is derived from the U.S.S.R. 
Cosmos satellite deceleration results collected over the 
period from 1964 through 1970 (Reference 21) , This model 
determines the atmospheric density by fitting density obser- 
vations over the range of altitudes and temperatures encoun- 
tered by the Cosmos satellites. Use of this model is 
restricted because the coefficients were empirically deter- 
mined over a limited altitude region and during only a por- 
tion of the 11-year solar cycle. The data were extended by 
using Jacchia's 1970 model (Reference 22). See References 
16 cuid 21 for further details. 

2.2.2 RECOMMENDED MODEL FOR THE UPPER ATMOSPHERE 

The Jacchia- Roberts model, with several modifications devel- 
oped in the Jacchia 1977 model, is recommended as the upper 
atmosphere model for LAIRS. The modifications involve using 
new parameters for the diurnal variation and the latitude- 
dependent time lag for the geomagnetic activity effect. With 
these modifications , the Jacchia-Roberts model explicitly 
includes most of the important atmospheric variations, and is 
not as complicated as the Jacchia 1977 or the MSIS model. The 
Jacchia-Roberts model can also be easily readjusted by updating 
only a few parameters which determine the basic structure of 
the temperature profiles, 

2.3 WIND MODELS 

Winds are a direct result of pressure and temperature differ- 
ences in the atmosphere and can be mathematically derived from 
a given pressure-temperature distribution. Alternatively, 
winds can be modeled by assuming a set of empirical wind pro- 
files as a default model and using observed wind velocities 
to adjust this model. The following sections elaborate on 
these two possibilities, drawing primarily on References 23 
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and 24 for the discussion of the analytic approach, and on 
the CIRA 1972 wind profiles (Reference 5) for the empirical 
model discussion. Figure 2-2 shows two typical wind pro- 
files. 

2.3.1 EVALUATION OF AVAILABLE MODELS 

2. 3. 1.1 Analytic Model 

The acceleration of an air parcel relative to the Earth can 
be expressed in terms of the pressure gradient and the accel 
eration due to the Coriolis force, frictional forces, and 
gravity. For frictionless flow, and neglecting the small 
terms introduced by the curvature of the Earth, the acceler- 
ation equations are as follows; 


^ = -ZcOuSlne - 4 - — 
dt f 


(2-29) 


(2-30) 


dtJ ^ L ^ _ a 

ddt P GZ J (2-31) 


where x, y, and z represent the eastward, northward, and 
upward directions, respectively, and u, v, and w are the. 
velocities corresponding to these directions. The angular 
velocity of the Earth is denoted by CO / 0 is the 

latitude, is the density, and p is the pressure. 
Because vertical wind velocities are much smaller than 
horizontal velocities, the Coriolis term has been deleted 
from Equation (2-31) . 
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Figure 2-2. Typical Summer and Winter Wind Profiles 
(30 Degrees North Latitude) 
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In the Earth's atmosphere, the accelerations on the left 
sides of Equations (2-29) through (2-31) are generally an 
order of magnitude smaller ' than the terms on the right . 
The geostrophic approximation consists of setting these 
accelerations to zero, so that the pressure gradient is 
balanced by the Coriolis force. This yields 

V =—L—iL 

ZfusinQ <dX 

u - 

Zf>u)sine ' 3 ^ 





(2-34) 


Equation (2-34) is the equation of hydrostatic equilibrium 
and Equations (2-32) and (2-33) are its analogs in the hori- 
zontal directions, i.e., the geostrophic wind equations. 

In the presence of a high or low pressure system. Equations 
(2-29) and (2-30) must be modified to include the centrifugal 
force term arising from air flow about the high or low. If 
R is the radius of curvature of the flow, the geostrophic 
approximation gives rise to the gradient wind, i.e.. 


V = "R sin 0 - J(-R cosine f- ^ |f 
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for flow about a high pressure system, and 


V = J(‘RcOsinef + -p- -Ro^S]r,9 (2-36) 

for flow about a low pressure system. In these equations, V 
denotes the magnitude of the horizontal wind velocity and 
Op/3n) denotes the radial gradient of the pressure. There 
is a limit to the velocity of flow about a high pressure 
system, since the term under the radical must be positive, 
but there is no corresponding limit to flow velocities 
about a low pressure system. The latter situation holds in 
the case of hurricanes . 

The geostrophic approximation can be used to derive 
(Reference 23) the vertical wind shear produced by nonuni- 
form horizontal heating: 


ivj (zcJslne) 



( 2 - 37 ) 


This is one form of the thermal wind equation. There is a 
similar equation for meridional flow (v) . Here R* is the 
universal gas constant and M is the mean molecular weight. 
Thus, given a horizontal temperature (T) distribution and 
the pressures p^ and P 2 at altitudes and (^2 ^ ^1^' 
one can determine the difference in wind velocity between 
those two heights. For most of the atmosphere, the combina- 
tion of the geostrophic wind (i^) and the thermal wind (i{^) 
has been found to be an adequate description of wind struc- 
ture, the thermal wind giving the increment in the geostro- 
phic wind with altitude, i.e., 

U^CZz^ - K-jC'Z,) -i- UrC^,,£z) 
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where tC-j* is the quantity given by Equation (2-37) and U.^ 
is given by Equation (2-33) . 

The geostrophic approximation, which gives rise to all of 
the above wind equations, is valid above the lowermost 1 or 
2 kilometers of the atmosphere (the planetary boundary layer) . 
There friction is important and topographical features make 
modeling of winds extremely difficult in general. Also, in 
the region between 80 and 120 kilometers altitude, ionospheric 
electro jet activity produces sporadic winds which cannot be 
treated by the wind equations (Reference 25) . Thus, the 
analytic model is of use only between the altitudes of 
(approximately) 2 and 80 kilometers. When the radius of 
curvatxire of the flow is known from sources such as satellite 
weather maps, the gradient wind equation can also be used. 

2 . 3 . 1 . 2 Empirical Model 

Vertical profiles of the zonal wind ( ti ) are given for 
various latitudes and for each month of the year by CIRA 
(1972) . These profiles cover the altitude range of 25-100 
kilometers . The meridional wind (v) , which is much smaller 
than the zonal wind at these altitudes, is not tabulated but 
may be estimated from accompanying diagrams. These profiles 
show significant monthly variations , mainly because of the 
strong jet stream around 60 kilometers which reverses direc- 
tion seasonally. This seasonal change of profile shape, 
together with the fact that above 80 kilometers the monthly 
profiles show little resemblance to each other, makes a 
polynomial fitting scheme impractical. Interpolation between 
data points stored on permanent files would have to serve as 
a default model. No model adjustment would be possible; the 
observed wind profiles would merely replace the default pro- 
files before the interpolation began. 
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The situation is much the same below 25 kilometers, except 
that there the source of the default data is information 
received directly from LaRC. As with the models for the 
other atmospheric parameters covered in Section 2.1, diurnal, 
semidiurnal, and latitudinal effects can be introduced into 
the wind model. The latitude variation can be derived from 
tabulated gradients of the wind components (based on infor- 
mation in References 5 and 24; the diurnal and semidiurnal 
variations can be modeled by Equation (2-15) . Amplitudes of 
the diurnal and semidiurnal east-west wind variations are 
charted in CIRA 1972 up to an altitude of 80 kilometers. 

Above that height, diurnal and semidiurnal wind variations 
have been studied by Fellous et (Reference 26) , who 

provide a mathematical scheme for approximating the east- 
west wind amplitude. It is known that the amplitudes and 
phases of the diurnal and semidiurnal wind variations change 
with season and may change from day to day, but these changes 
have not been well studied. 
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2.3.2 RECOMMENDATION OF A WIND MODEL 

As pointed out in Section 2. 3. 1.1, the analytic model of wind 
structure is limited to regions below an altitude of 80 
kilometers. The primary advantage of the analytic model 
is that it requires no input other than temperature, pressure, 
and density distributions at a given time and location. If 
longitudinal, latitudinal, diurnal, and seasonal effects have 
been adequately modeled in the state variables, they will 
automatically be included in the calculated winds. However, 
enough measurement profiles to provide a good estimation of 
the latitudinal and longitudinal (or hour angle) gradients 
of temperature and pressure are required by this model. It 
will not be adequate to use default values of these gradients 
together with measured profiles to obtain the wind components , 
as these components are linearly proportional to the gradients ; 
i.e., a 100-percent error in a pressure gradient would produce 
a 100-percent error in calculated wind speed, although it 
would produce only a 1-percent or 2 -percent error when used to 
extrapolate pressure a few degrees (most such gradients are 
quite small) . In addition, sporadic gusts and turbulence 
cannot be modeled by the wind equations . 

The empirical model has important drawbacks as well. One is 
that no actual modeling of vertical wind structure is possible — 
only interpolation between default or measured data points . 
Another drawback is the incomplete knowledge of diurnal and 
semidiurnal wind variatons. These must be included not only 
because the Shuttle covers several hours of local time in a 
few minutes , but also because there may be a time lag between 
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the measurement of the wind profiles and the reentry of the 
Shuttle (because daily variations are quite small below an 
altitude of 30 kilometers and because the Shuttle spends 
little time and covers little horizontal distance in this 
.region, it is recommended that diurnal and semidiurnal varia- 
tions not be modeled below 30 kilometers) . 

In summary, either the analytic or the empirical model can 
be used to model the wind. If there will not be enough 
meteorological data available to allow accurate determination 
of temperature and pressure gradients, the analytic model 
should be avoided. Therefore, it is recommended that the 
empirical model serve as the primary method of wind deter- 
mination, since it is not yet certain how many meteorological 
profiles will be available for a typical Shuttle reentry. 
However, certain routines should be provided so that the 
analytic model can be used in the event that accurate deter- 
mination of pressure and temperature gradients becomes pos- 
sible . 
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SECTION 3 - FITTING PROCEDURES 


Meteorological data obtained from the National Oceanic and 
Atmospheric Administration (NOAA) contain vertical profiles 
of temperature/ pressure, density, and wind at a specified 
longitude, latitude, and local solar time. This section 
presents a description of a method for finding adjusted 
model altitude profiles, adjusted diurnal and semidiurnal 
variations, and the adjusted latitude dependence of the atmos 
spheric model from these meteorological data. Section 3.1 
contains a discussion of the diurnal and latitudinal varia- 
tions and Section 3.2 deals with the fitting of altitude 
profiles. Section 3.3 dicsusses wind profiles and 
Section 3.4 contains some discussion of accuracy estimation 
for the fitting procedures. 

3.1 DIURNAL, SEMIDIURtlAL , AND LATITUDINAL VARIATIONS 

The observed vertical profiles will be labeled by the local 
solar time, t, and the latitude, 0. Local solar time 
effects, i.e., the diurnal and semidiurnal variations, and 
the latitude dependence of the atmosphere will be explicitly 
considered insofar as the distribution of observation data 
allows. New diurnal and semidiurnal amplitude and phase 
relationships can be obtained from a set of altitude pro- 
files, all having the same reference latitude, 9„. These 

r\ 

profiles can be obtained by latitude translation (to latitude 
using the default (CIRA) latitude dependence. These 
profiles must have a reasonably good distribution in local 
solar time so that new diurnal and semidiurnal variations 
can be calculated. Having found a new diurnal and semi- 
diurnal variation, a new latitude dependence can be deter- 
mined from a similar set of vertical profiles, all having 
the same reference local solar time, tj^, with a finite 
spread in latitude. This set of vertical profiles can be 
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constructed by local solar time translation using the new 
diurnal and semidiurnal variations. Finally, a large 
composite reference altitude profile at ( tj^) can be 

constructed from the original set of observed vertical pro- 
files using the newly obtained diurnal and latitudinal varia- 
tions . 

_ • 

Let ^ ) denote the ith vertical profile in the 

original set of observed vertical profiles, and let 

and L. ( ©* , 0/^ ) denote the newly computed diurnal and 

latitudinal translations, respectively. First 'P^ ( 9^ 

is translated from ( P. , t.) to the reference point 
/\ 1 i • 

( • This translated profile is denoted as 1^*" ( &ji 

« • 

Then the reference profile, , is constructed as 


( 3 - 2 ) 

i. 

This reference profile is a composite set of many profiles 
merged into one. 'PC&R ,t.R ) will be used to fit a model 
altitude profile at ( , tT/? ^ 

3.2 FITTING PROCEDURES FOR VERTICAL PROFILES 

Reference vertical profiles constructed in Section 3.1 for 
temperature, pressure, and density will be used to derive 
respective model altitude profiles. Temperatures will be 
fitted to a set of polynomials of Z for the lower atmosphere 
and to the Jacchia-Roberts temperature profile for the upper 
atmosphere. In the case of pressure and density, C P/ pc') 

and &n ( P/a ) will similarly be fitted to a set of 
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polynomials and to the Jacchia-ROberts model. The quantities 
p and P denote the pressure and density^ respectively, and 
the subscripted quantities represent values at a reference 
altitude. 

3.2.1 TEMPERATURE PROFILES 


The entire altitude range is divided into four regions as 

shown in Figure 3-1. At the boundary between two consecutive 

regions , the temperature is required to be continuous . The 

Jacchia-Roberts temperature profile will be used for the 

B 

uppermost region which starts at an altitude (Z^) near 120 
kilometers. When there are enough observation data available 
to cover up to and beyond 125 kilometers, will be taken as 
125 kilometers. Then the Jacchia-Roberts temperature profile 
to be used for the uppermost region (Z^125 kilometers) takes 
the following form: 


Ti.(z) = %-(T^-T,^)e 


-BCZ) 


(3-3) 


H(Z) = (?® 2. -Zf ^ Zf 


where 

t; 

77 ® 


exospheric temperature 
temperature at Z® 

a 

upper temperature gradient at z” 
constant = 6356.766 kilometers 
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Figure 3-1. Temperature-Altitude Regions 
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For altitudes below 125 kilometers, the following sets of 
polynomials will be used: 


T-Cz) Z c,„fz-zf„y 


^ z. yz 


( i= Z, 3, 


(3-4) 


where 

_ — B 8 

= temperature at 

= order of the polynomial for region i 

Parameters * 7 ^^ , C^n 5 i = 2., 3, 4^ } 

and a part of can be readjusted such that Tj_(Z) (i=l'^4) 

best fit the observed temperature data. When there are not 
enough data to cover all regions, a subset of the above param- 
eters can be selected to be solved for. Mathematical details 
of these procedures will be given in Section 4.3.2. 

3.2.2 DENSITY AND PRESSURE 

Pressures and densities could be obtained from the tempera- 
ture profile computed in Section 3.2.1 using Equations (2-13) 
and (2-14) . However, the quantities can also be fitted 
directly if observation data for temperature, density, and 
pressure are available from independent sources. In this 
case, the entire altitude range is divided into four regions 
as in the case of temperature. The Jacchia-Roberts density 
model will be used for the uppermost region and a set of 
polynomials for the lower regions. A set of functions 
{F^(Z)} defined as 
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will be used in the fitting procedure. The ^(Z) term 
denotes a density value at Z and denotes a reference 

density at a reference altitude. 


Fi(z) =-irt ) Z5zf 

Rc(Z)= F.® i-lcnCZ-Zi,)” 

r\=f 

zf., > z > z® 

( i = Z, 3, ^ ) 

where 

Jacchia-Roberts density 


Parameters F? (i = 1,2,3), C. (n=l~N.; i=2/v4) and the 
1 ' ' ' in 1 ' 

Jacchia-Roberts temperature parameters can be adjusted such 
that best fits the observed altitude profile for 

the atmospheric density. 


A similar approach can be used for pressure. In this case 
Fj^(Z) is defined as follows. 


Fjz) = -U 
P„(Z)= 

^ ^ ^ : Z. y Ido Km. 
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In these equations, M(Z) denotes the mean molecular weight, 
the number density of each individual species, and 
Avogadro's niimber. The first expression for Pj^^ applies to 
the region of barometric equilibrium and the second expression 
applies to the region of diffusive equilibrium. 
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3 . 3 WIND 


Observed wind data will be processed according to the same 
procedures described in Section 3.1 for the diurnal, semi- 
diurnal, and latitudinal variations of temperature, pressure, 
and density. However, its vertical profile will not be 
fitted to any model, but rather a simple and straightforward 
interpolation scheme will be used throughout to find wind 
velocities at various points along the Shuttle trajectory. 

As an option, a theoretical wind profile can be created 
using the spatial variations of temperatures and pressures 
(see Section 2.3). 

3.4 ACCURACY ESTIMATION 

Examples of the NOAA meteorological profiles that are expected 
to be available after each Shuttle flight have been provided 
by LaRC.^ The measurement errors accompanying these profiles 
show that, in the lowest 25 kilometers of the atmosphere, the 
uncertainties in temperature, pressure, and density are all 
less than 2 percent. Between 25 and 75 kilometers, they may 
reach 10 percent or more. Numerical experiments using this 
sample meteorological data and the polynomial fitting proce- 
dure discussed in Section 3.2 indicate that errors resulting 
from the polynomial fitting itself can easily be kept within 
this range. For example, fitting two third-order polynomials 
(constrained to join and be continuous at an altitude of 
25 kilometers) to the sample data between the surface and 
65 kilometers results in a standard deviation of 4° K, or 
about 1.5 percent of the typical temperature in that region. 
Even better results (1° K - 3° K) are obtained when the 
smoothed CIRA or ERA profiles are fitted. 


1 


Private communication from J. M. 


Price . 
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The situation with wind errors is less encouraging. The 
sample meteorological profiles show that errors of between 
10 and 20 percent of the total wind speed are possible below 
75 kilometers, and that uncertainties of up to 100 percent 
may occur above that level. No uncertainties or standard 
deviations are given for the CIRA. profiles. However, mean 
deviations of observed winds from modeled winds are given 
for certain locations. These deviations are quite variable, 
but generally show the same situation as do the meteorological 
data: below 75 : kilometers , the modeled winds are often off by 

10 or 20 percent; above 75 kilometers they are sometimes off 
by as much as 100 percent. These errors must be accepted as 
the limits of wind modeling. When wind uncertainties are 
considered, however, it should be recalled that in. the alti- 
tude region where the largest errors are expected, the Shuttle 
velocity will be greater than 5 kilometers per second. The 
wind speeds there are at most 100 meters per second, so even 
if the winds were 100 percent in error, this error would be 
no more than 2 percent of the Shuttle velocity. 

It is also important to consider errors in the modeling of 
diurnal and semidiurnal variations. As pointed out previously, 
below 25 kilometers, these effects are small enough to be 
ignored, except in the planetary boundary layer. Above 25 
kilometers they grow rapidly, primarily because of absorp- 
tion and re-radiation by the ozone layer. The amplitude of 
the diurnal temperature tide can reach 10° K (4 percent) 
around 50 kilometers; the amplitudes of the diurnal and 
semidiurnal components of pressure and density variations 
range between 6 percent and 12 percent near 80 kilometers. 

Wind amplitudes at such heights may reach 10 to 30 meters 
per second. Thus, in 3 to 6 hours (local time) the atmos- 
pheric parameters undergo variations that can exceed measure- 
ment error at some altitudes , and should therefore be 
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incorporated into the model. Unfortunately, the amplitudes 
and phases of these variations change greatly with the 
seasons, and these changes have been poorly studied. 

Sources such as CIRA 1972 provide only an incomplete mapping 
of diurnal and semidiurnal effects. Reliance on the default 
model of daily variations (as would be necessary when only 
one or two profiles were available) could conceivably intro- 
duce errors of 10 to 12 percent in some parameters above 50 
kilometers. 

The effect of errors in the latitude gradients may be readily 
estimated. From the CIRA 1972 profiles it can be seen that 
these quantities are quite small, with typical temperature 
gradients at 30 degrees north latitude being only a few 
tenths of a degree Kelvin per degree latitude. Pressure 
and density gradients are marginally more significant. Still, 
it is extremely unlikely that errors in the values of the 
latitude gradients would produce uncertainties in the model 
parameters exceeding the uncertainties due to measurement 
error. 

Finally, the consistency of the models studied will be briefly 
considered. As discussed in Section 2.1, the difference 
between the average of the twelve 30-degree-north-latitude 
CIRA temperature profiles and the mean ERA 1975 profile is 
less than 1° K at 25 kilometers. The differences for pressure 
and density are greater, but are still less than 4 percent. 
These are about the same magnitudes as the differences 
between either of these models and USSA 1962 near 25 kilom- 
eters. The Jacchia 1977 model assumes a temperature mini- 
mum at 90 kilometers, and takes the USSA 1976 value of 188° K 
at that point. The CIRA profiles at 30 degrees north aver- 
age 189° K at 90 kilometers. The difference between the CIRA 
density and the Jacchia 1977 (USSA 1976) density at 90 
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kilometers is 3 percent. Thus, all the models considered 
are consistent. The combination of the ERA 1975, CIRA 1972, 
and Jacchia-Roberts models has been recommended because it 
is believed to offer the most comprehensive, flexible, and 
easily implemented atmospheric modeling scheme. 
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SECTION 4 


FUNCTIONAL SPECIFICATIONS 


This section details the software requirements and specifica- 
tions for the Langley Atmospheric Information Retrieval Sys- 
tem (LAIRS) . Mathematical specifications are also discussed, 
along with input/output requirements . 

4.1 SOFTWARE REQUIREMENTS 

The LAIRS software must be designed to perform the following 
functions ; 

• Respond to user input specifying program options 

• Receive an input ENTREE Program trajectory file 
containing Shuttle position versus time 

• Receive input meteorological measurements contain- 
ing vertical profiles of temperature, pressure, 
density, horizontal wind speed, and horizontal wind 
direction from various stations along or near the 
Shuttle reentry path 

• Produce values of the atmospheric parameters (temp- 
erature,. ^ pressure, density, and horizontal wind com- 
ponents) for points on the Shuttle trajectory, 
using the input meteorological data as a basis for 
calculation 

• Provide default values of the atmospheric parameters 
when no or insufficient meteorological observations 
are available 

• Produce output reports and an output user file con- 
taining the values of the calculated parameters 
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The calculation of the atmospheric parameters can take one 
of three forms: calculation based on the default coefficients 

derived from the Jacchia-Roberts , CIRA 1972, and ERA 1975 models 
discussed in Section 2; calculation based on adjusted coeffi- 
cients determined from the default coefficients and the input 
meteorological data by the fitting procedure discussed in 
Section 3; and a simple interpolation of the nearest meteoro- 
logical profile with no attempt at any sort of modeling. 

The default coefficients will be stored on permanent files 
and the adjusted coefficients will be stored on working files. 
The software must thus be able to create and access the 
appropriate files in order to perform the necessary calcula- 
tions . 

The software for adjusting the model coefficients will be 
the most complex part of the system. As discussed in 
Section 3, least-squares techniques will be the basis of 
the adjustment procedure. It will be possible, however, 
to adapt certain ENTREE Program routines to the model-fitting 
process, thereby lessening the development time required for 
this section of LAIRS. 

The following options will be available to the user: 

• Calculate parameters based on the default model, 
on the adjusted model, or on the interpolation 
scheme 

• Produce an Atmospheric Parameter Output Report 
(APOR) giving the calculated parameters along the 
Shuttle trajectory 

• Produce a Local Atmospheric Profile (LAP) version 
of the Atmospheric Parameter Output Report, giving 
a set of vertical profiles of the calculated param- 
eters for a specified latitude, longitude, and time 
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• Rebuild the permanent files using input meteorologi- 
cal observations . This option would be useful in 
the event new CIRA profiles become available, for 
example . 

4.2 INPUT 

4.2.1 TRAJECTORY DATA 

The LAIRS input trajectory data will be the output trajectory 
file of LaRC's ENTREE Program, giving Shuttle position as a 
function of time. The contents of this file that will be 
utilized by LAIRS are shown in Table 4-1. 

4.2.2 OBSERVATION DATA 

The meteorological observations required by the model adjust- 
ment process will be supplied by NASA on a nine-track magnetic 
tape after each Shuttle flight. The format will be that of 
the Meteorological Input Tape described in Reference 21 . The 
information supplied will consist of a set of vertical pro- 
files of temperature, pressure, density, wind speed, wind 
direction, and the uncertainty in each of these measurements 
(the dewpoint is also provided on the tape, but will not be 
used by LAIRS) . These profiles will cover the region from 
0 kilometers to at least 120 kilometers, and are expected 
to be provided by several stations along the Shuttle reentry 
path. Each profile is assumed to have a time tag (GMT) that 
is within one-half hour of the time the profile was obtained. 
It will not be necessary for each point in the profile to be 
tagged. The foimtiat of the meteorological data required by 
LAIRS is shown in Table 4-2. ^ 

4.2.3 CONTROL DECKS 

LAIRS will be designed to accept as input a limited number 
of keyword cards giving the specific options selected from 
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Table 4-1. 


Variables to Be Read From ENTREE Trajectory File 


VARIABLE 

TYPE 

LENGTH 

DESCRIPTION 

TIME 

REAL 

1 

Time tag of state vector 
(GMT) 

X 

REAL 

3 

State vector: altitude 

(meters) , latitude 
(degrees east) , longi- 
tude (degrees north) 

XSTD 

REAL 

3 

Standard deviations of the 
state vector components 
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Table 4-2. 


Format of LAIRS Meteorological Data Tape 


(1 of 2)^ 


WORD 

SYMBOL 

DESCRIPTION 

UNITS 

1 

LAT 

Latitude 

Degrees, + N 

2 

LON 

Longitude 

Degrees, + E 
to 360 

3 

FLAG 

Data flag; 

= 0, measured data 
= 1, modeled data 
= 2, combined measured 
and modeled data 


4 

— 

Spare 


5 

ALT 

Geometric altitude 

Feet 

6 

WS 

Horizontal wind speed 

Feet/second 

7 

WD 

Direction horizontal 
wind is coming 
relative to true 
north, north being 
0 degrees, increas- 
ing positively 
clockwise 

Degrees 

8 

TE 

Ambient temperature 

Degrees C 

9 

PR 

Ambient pressure 

Millibars 

10 

D 

Ambient density 

Grams/meter 

11 

DW 

Dew point (not used) 

Degrees C 

12 

TEU 

Ambient temperature 
systematic uncer- 
tainty 

Degrees C 

13 

PRU 

Ambient pressure 
systematic uncer- 
tainty 

Millibars 

14 

DU 

Ambient density 
systematic uncer- 
tainty 

Grams/me ter^ 

15 

HWSUS 

Horizontal wind speed 
systematic uncer- 
tainty 

Feet/second 

16 

HWSUN 

Horizontal wind speed 
noise or fluctua- 
tion uncertainty 

Feet/second 


^All quantities are real variables one word in length. 
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Table 4-2. Format of LAIRS Meteorological Data Tape (2 of 2) 


— 

WORD 

SYMBOL 

DESCRIPTION 

UNITS 

17 

VWSUN 

Vertical wind speed 
noise or fluctua- 
tion uncertainty 

Feet/second 

18 

HWDUS 

Horizontal wind direc- 
tion systematic 
uncertainty 

Degrees 

19 

HWDUN 

Horizontal wind direc- 
tion noise or fluc- 
tuation uncertainty 

Degrees 

20 


Spare 

— 
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those listed in Section 4.1. The information on these cards 
will include the following: 

• Model to be used (default, fitted, or interpolation) 

• Type of calculation profile (ENTREE Program trajec- 
tory, or Local Atmospheric Profile, or both) 

• Decision to rebuild permanent files 

Other possible options and the models to which they apply 
are summarized in Table 4.3. 

Also given on keyword cards will be specifics related to 
various options (such as LAP location) . A set of station 
cards specifying latitudes and longitudes of stations 
supplying meteorological profiles, together with the number 
of profiles supplied by each station, will be required. 
Finally, the option to rebuild the permanent files will 
require special input similar to the meteorological obser- 
vations (the exact nature of this input will depend upon 
the specific routines devised for model fitting) . 

4.2.4 PERMANENT AND WORKING FILES ; ' 

As mentioned in Section 4.1, it will be necessary to main- 
tain both permanent files, which are associated with default 
model parameters, and working files, which are associated 
with adjusted model parameters, for the operation of LAIRS. 
The permanent files will be accessed by the model-fitting 
routines only. These routines will create the working files, 
which will be accessed by the parameter calculation routines. 
When the default model is requested, the model-fitting rou- 
tines will take the appropriate information from the perma- 
nent files and write it to the working files unaltered. 

The content of the files will be described in the following 
sections. 
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Table 4-3. Options Related to Model Choice 


— MODEL 

OPTION — 

DEFAULT 

FITTED 

INTERPOLATED 

Altitude region in which to use model 
(default = entire) 


/ 

/ 

Number of default data points to be added 
between Z1 and Z2 (default = none) 


/ 


T, P, p fitting done independently 
(default) or P and p derived from T 


/ 


Method to be used in joining Jacchia- 
Roberts to lower atmosphere 

/ 



Inclusion of wind (default = yes - 
empirical) 

/ 

/ 

/ 

Use of analytic wind 


v/ 


Use of diurnal and semidiurnal variations 
and latitude gradients 

/ 

/ 

/ 




4. 2. 4.1 PERMANENT FILES 

For the lower atmosphere, the permanent files will consist 
of several sets of polynomial coefficients representing the 
vertical distribution of temperature, pressure, and density 
for each month of the year. Horizontal wind components and 
latitude gradients of the atmospheric parameters will be 
tabulated on the files at 5-kilometer intervals for each 
month. Diurnal and semidiurnal variation coefficients will 
also be provided, but they will be identical for all months. 


The Jacchia-Roberts model recommended for the upper atmos- 
phere will require three coefficients (see Section 3.2). 
These will not have to be tabulated for each month. Addi- 
tionally, the Jacchia-Roberts model requires three other 
files that are not, strictly speaking, permanent, but are 
also not derived from the meteorological observation files. 

These are the 3-hour geomagnetic planetary index (K ) file, 



the 10.7-centimeter solar flux (Fj_q - j ) file, and a file 
consisting of an 81-day running average of 7 ^ ' 

centered on the day of interest. These files must have 
new values added periodically, but once new values are 
added, they remain on the files permanently. 


4. 2. 4. 2 WORKING FILES 


The working files for the upper and lower atmospheric models 
will contain exactly -the same parameters as the permanent 
files, except that the working files will be produced by the 
model-adjustment software using the meteorological input 
data and the fitting processes discussed in Section 3. There 
will be no working files corresponding to the Jacchia-Roberts 

^p' ^10.7' ^10.7 
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The data contained in the working files are summarized in 
Table 4-4. The content and format of the permanent files 
will be identical to that of the working files, except that 
there will be 12 sets of permanent files (one for each month) . 
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Table 4-4. Contents of LAIRS Permanent and Working Files 


LOWER ATMOSPHERE: 


Three sets of polynomial coefficients for each of T, 
P, and p (0-25, 25-65, and 65-110 kilometers) 

Tabulated winds U and V (0-110 kilometers) 

Tabulated diurnal and semidiurnal coefficients for 
each of T, P, p, U, and V (0-110 kilometers) 


Tabulated latitude gradients for each of T, P, p, U, 
and V (0-110 kilometers) 

UPPER ATMOSPHEPvE: 

2 B B 

Jacchia-Roberts coefficients T , T, , and G. 

oo' 1' 1 

Tabulated wind U and V (110-130 kilometers) 


Geomagnetic index^ 

3 

Solar flux - F^q ^ 

3 

Average solar flux 



- F 


10.7 


^Permanent files will consist of 12 monthly sets of default 
values for the quantities tabulated here. 

2 

The Jacchia-Roberts model includes diurnal, latitude, and 
annual effects. Only one set of default values will be 
required. 

3 _ 

Kp, Fjj^q and F^^ ^ need to exist only on permanent files. 
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4.3 MATHEMATICAL SPECIFICATIONS 


Mathematical details necessary for the LAIRS implementation 
will be presented in this section. The description of the 
default atmospheric model is given in Section 4.3.1, followed 
by the discussion of fitting procedures in Section 4.3.2. 
Section 4.3.3 contains a summary of the computational proce- 
dures for various atmospheric parameters. 

4.3.1 DEFAULT ATMOSPHERIC MODEL 


The default atmospheric model for LAIRS consists of three 
parts: the Edwards Air Force Base Reference Atmosphere from 

the ground to approximately 25 kilometers, the Cospar Inter- 
national Reference Atmosphere (CIRA 1972) from 25 kilometers 
to 110 kilometers, and the Jacchia -Roberts atmospheric model 
above 120 kilometers.* 


The Jacchia-Roberts atmospheric model determines the atmos- 
pheric mass density, p, and the number density of individual 
species, n^, using the following two equations: 



: : Z ^ 100 Km 


(4-1) 


dz. 




Me 3 

K*T 


Z >100 Km 


(4-2) 


where = thermal diffusion coefficient of the species i 

Mj^ = molecular weight of the species i 
M = mean molecular weight 
g = gravitational acceleration 
R* = universal gas constant 


*An interpolation scheme will be used between 110 and 120 
kilometers. 
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The first equation given above is the barometric equation 
and applies to a region of the homogeneously mixed atmosphere 
below 100 kilometers and the second equation is the diffusion 
equation for the number density of each individual species 
and applies to the atmosphere above approximately 100 kilo- 
meters, where the lighter species separate from the heavier 
species by the process of molecular diffusion. 

The temperature, T, in Equations (4-1) and (4-2) has to be 
known as a function of Z. Temperature-height profiles were 
empirically obtained in the Jacchia 1971 model and will be 
slightly modified according to the formulation given by 
Roberts so that the density equations can be integrated an- 
alytically (References 15 and 18) . 


The temperature profile to be used in the Jacchia-Roberts 
atmospheric density model takes the form of a fourth-order ' 
polynomial for heights from 90 kilometers to 125 kilometers 
and an exponential function of Z for heights above 125 kilo- 
meters 


T(z) = T, f(r^-T.)Z.a„vz.-Zxf 


n=i 


(4-3) 


This equation gives temperatures between Zq and Z^, and the 
coefficients a^ are determined from the following conditions; 


T(Zq) “ "^0 ■ temperature 

T(Z ) = T : temperature at inflection point 




= a^(T^-Tj) 


(4-4) 


1.9 


(Tx-Tq) 



an empirical 
relation 
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0 


(4-4) 
(cont ' d) 



2a2(T^-To> = 


The four coefficients (i = l-»-4 ) can be expressed in terms 
of Zq, Z^, Tq, and as follows: 


a. = /.9 /(Z>,-z,) 

dz - 0 

Uj =-/.7/('z,-Zj* 
cif. =-«-sAz,-zj« 


The values of Zq and Z^ are fixed at 90 kilometers and 125 

kilometers, respectively. The value of Tq was 183° K in 

the Jacchia 1971 model and was changed to 188° K in the 

Jacchia 1977 model. The quantity T is determined from the 

value of T (Reference 13) 

00 ' ' 



^ -f 


^71 

CD 


•f c e 


IclZ 


Ob 


where a 
b 
c 
k 


371.6678 

0.0518806 

-294.3505 

-0.00216222 


(4-6) 
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For heights above 125 kilometers, Roberts used the following 
function to define the temperature: 


TCZ) = - (Ta-Ty) 



Ra.'^Zi / 


(4-7) 


where G = temperature gradient at Z 

X X 

R, = 6356.766 kilometers = mean radius of the Earth 

A temperature profile defined by Equations (4-3) and (4-7) 
is continuous at Z = Z regardless of the choice of G . The 
temperature gradient will be continuous at Z = Z if 

X 

G = 1 . 9 (T -T^) / ( Z -Z» ) . However, in the Jacchia-Roberts 

X X U X u 

density model implemented in GTDS, G„ was determined as a 

X 

function of such that the resulting density profiles 
obtained using Roberts' temperature profiles gave the best 
least-squares fit to Jacchia's tabulated density data. In 
this case, the temperature gradient is not continuous at 
Z = Z . 

X 

In addition to temperature profiles. Equations (4-1) and 
(4-2) require a knowledge of g(z) and "mCz) . i The following 
expressions were used for these quantities:; 


2 

where gQ = 9.80665 meters per second , and 

— ^ n 

M CZ1) = 2_4 Cn ( Z- •• 9o<z<.ioo (4-9) 

n=o 

C Z In fCm ) 
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Numerical values in this equation are given in Reference 14. 
Roberts analytically integrated Equations (4-1) and (4-2) 
using temperature profiles given by Equations (4-3) and (4-7) ; 
the gravitational acceleration, g(z); and the mean molecular 
weight, M(Z) . The quantities g(Z) and M(Z) are defined by 
Equations (4-8) and (4-9) , respectively. The solutions ob- 
tained by Roberts, however, contain rather lengthy algebraic 
expressions and will not be given in this report. (See Ref- 
erences 15 and 18 for further details.) 

The computation of T^ is described in Section 2.2. Discus- 
sion of various other corrections included in the Jacchia- 
Roberts density model can be found in References 13 and 15 . 

The Jacchia-Roberts atmospheric density model will be used 
for heights above 120 kilometers and the CIRA model for 
heights below 110 kilometers. For heights between 110_ kilo- 
meters and 120 kilometers, an appropriate interpolation 
scheme will be used for each atmospheric parameter. A four- 
point Langrange interpolation method will be used for temp- 
eratures and a two-point exponential interpolation scheme 
for densities and pressures. The latter interpolation pro- 
cedure is defined as 

(Z-Z,)] (4-10, 

which gives p = at Z = Z^ and p = Pg at Z = Zg. 

An alternative method of connecting the CIRA model to the 
Jacchia-Roberts model involves readjusting some of the pa- 
rameters that determine the structure of the Jaccia-Roberts 
model. Let Z^ be the boundary between the two regions de- 
scribed by the CIRA and Jacchia-Roberts models, respectively. 

Z- can assume any value between Z (90 kilometers) and Z 
1 o X 
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(125 kilometers) , excluding Z . Two parameters, T and p , 

X o o 

in the Jacchia-Roberts model can be readjusted so that the 
various atmospheric parameters will be continuous at Z = Z^. 
Let and be the newly adjusted values of these two 
parameters. Then, using Equation (4-3) , 

X 

n-\ 

For calculating the density parameter, the following simple 
scaling law will be sufficient. 

fo' tc Pj^CZ,) 

where i^(Z,) should be calculated using T^. 

JR 1 o 

When the temperature and density are continuous at Z = Z^, 
the pressure will also be continuous, at least approximately. 
In this alternative method, there is no need for any inter- 
polation scheme. The region above Z = Z^ is described by 
the Jacchia-Roberts model and the region below Z = Z^ is de- 
scribed by the CIRA model. 

If Z^ is chosen to be near 90 kilometers, then it is more 
consistent and advantageous to completely redetermine Z^, 

T^, and p^, i.e., the position, temperature, and density of 
the last temperature minimum, from the default CIRA tempera- 
ture profile. In this case, will be used in place of Z^. 

The CIRA 1972 model will be used in LAIRS to compute various 
atmospheric parameters for heights between 25 and 110 kilo- 
meters, and the ERA model will cover the region from the 
ground up to 25 kilometers. Temperatures, densities, and 
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pressures near 30 degrees north latitude will be computed 
from a set of polynomials defined for each month of the year. 

The form of these polynomials was defined in Equations (3-4) 
and (3-6) in Sections 3.2.1 and 3.2.2, and the corresponding 
coefficients will be stored in the default atmospheric file. 

In addition to these coefficients, the default atmospheric 
file will also contain latitude gradients of various atmos- 
pheric parameters near 30 degrees north latitude and will 
contain amplitudes and phases of the diurnal and semidiurnal 
variations . 

4.3.2 DEFINITION OF FITTING PROCEDURES 

4. 3. 2.1 Diurnal, Semidiurnal, and Latitudinal Variations 

A new set of amplitudes and phases of the diurnal and semi- 
diurnal variations can be determined if observed data have 
enough spread in local solar time. Similarly, a new set of 
latitude gradients of atmospheric parameters can be deter- 
mined if the observed data have enough spread in latitude. 

The ith vertical profile in the original meteorological 

• » • . 

measurements will be denoted as 

The label i belongs to a set Jq which contains all labels 
of the profiles available in the original meteorological 
data, and Mq is the total number of meteorological measure- 
ments included in the ith profile. The quantity ? is used 
to denote an atmospheric parameter (temperature, pressure, 
or density) . 

The first step in finding new diurnal and semidiurnal vari- 
ations involves computing C (Z„,0i,,t.) from ? (Z ,0.,t.), 
using the default latitude gradients of ? and an appropriate 
interpolation scheme for Z. The parameters 0^^ and {Z^^} denote 
a reference latitude and a set of reference vertical points, 
respectively. 
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If •[4^(Zm,0c}tiy I j denotes the interpolated 

profile, then the latitude translation is performed on this 

i R 

interpolated profile to obtain S (Z^,6„,t.). This latitude 

in K 1 

translation is given as 




where g(Z ,6.,t.) denotes the latitude gradient of C at 

In the latitude translation, the second-order and higher 
order terms in (6 td-0i) and the local solar time dependence 
of g will be neglected. Within this approximation, 9 
can be replaced by 9 q(Zj^) , the default gradient at some ref- 
erence latitude. Interpolating ? from {Z^} to {Z^}, a three- 

m m 

point Lagrange method will be sufficient. 

The interpolation scheme should be used only in the neighbor- 
hood of real observation points. The number of measurements 
R 

at Z with different local solar times may not be a constant, 
^ R 

but may depend upon Z^. That is, at different heights, ob- 
servation data might have different distributions in local 
solar time.- This does not present any real problem as far 

as determining new diurnal and semidiurnal variations, since 

R 

these variations are determined at each height Z^ independ- 
ently. 

Using this set of profiles defined at 0^^ and a new set 

of diurnal and semidiurnal amplitudes and phases will be 

determined as follows. First, ^(Z^,6„,t) is assumed to have 

m 

the following t-dependence : 

ZS, ©K.t ) = /l.{ 1 1 - R, COifut TyS, ) +/? 2 cos Czwt } (4-12) 
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where to is the rotational speed of the Earth. The diurnally 
averaged amplitude, A^; the diurnal parameters, and 3^; 

and the semidiurnal parameters, R~ and 3_, are functions of 

R ^ ^ 

Z and only and can be determined using the set of profiles 

LJ 

constructed at and { . The set of equations that deter- 
mines these parameters is given by 


t.- ) = /l, f 1 f /?, WJ t 1- RjCoS 3 ) 

( ie j) 

As noted earlier, the number of profiles, N^., contained in 

R ^ 

J could be a function of Z^. If Nj>>5, the five parameters 
Aq , R 2 , 3j , S >2 can be obtained by least-squares fitting 

procedures. If Nj = 5, these parameters can be algebraically 
solved for. If Nj < 5, only a subset of these parameters can 
be determined. 

However, Nj is not the only factor that determines the solv- 
ability of Equation (4-13) . There should also be enough 
spread in local solar time. If all the profiles belonging to 
J are clustered within a couple of hours, not all of these 
parameters can be accurately solved for. In this case, a 
partial readjustment of these parameters will be carried out. 

After these diurnal and semidiurnal parameters are determined, 
new latitude gradients will be estimated. This involves com- 
puting 5^(Z^,0.,t„) from (Z^, 0 . ,t . ) , using the newly deter- 
mined diurnal and semidiurnal parameters as follows; 




(4-14) 
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(4-15) 




I -f ??1 cosCcotf? t/5,) + (z cotfi t^z) 

i -t R, cos(u)ti -fyS,)+ KiCos 


i R 

If the set of profiles, ^ enough spread in 

latitude in the neighborhood of the reference latitude, 6^^, 
the 6-dependence of ^ can be determined from 




The 6-dependence of C is defined as 


The two quantities 5(Z^,6_,t„) and g(Z^,6„,t„) will be deter- 

in x\ x\ m i\ i\ 

mined from the following set of equations; 

|'(z^,0^,tR)=§fzle,t^)+j(zS,e/!,t,:)K--eR) (4-17) 

If the nxamber of points in J, Nj, is greater than 2, least-| 
squares fitting procedures can be used to find ^(Z^,6_,t„) 
and g(Z^,6„,t„). If = 2, these two quantities can be 

III K K u 

obtained algebraically. 


The composite reference profile defined in Equation (3-2) 
can then be built from the original set of data as follows: 



i‘'(zi,,&x,u) = %^(zi,ei.u)tjcK,d,-,t,)(e^-d,) (4-i9> 


Here D(t.,t_) and g(Z^,0.,t.) are the new diurnal, semidiurnal, 
IK mil 

and latitudinal variations determined above. It should be 
noted that D and g are defined at . An appropriate in- 

terpolation scheme can be used to compute D and g at 

The composite reference profile to be used in the vertical 
fitting procedures takes the following form: 


The number of observation points in this super set is 

M = Z 

t = i 

Using M, the entire set can be written as 


^ (- ) I m-= 1''' m"! 
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4. 3. 2. 2 Vertical Profiles 

Model vertical profiles for temperature, density, and pres- 
sure use the Jacchia-Roberts density model in the uppermost 
region above 125 kilometers and a set of polynomials for 
heights below 125 kilometers. To define these functions, it 
is convenient to use a step function defined as follows : 


h(?) = 1 if ?>0 
= 0 if ?<0 


(4-21a) 
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(4-21b) 


Hi.(Z) 


h(Z-Z^) h(Zj-Z) 
1 if Z^ < Z < Zj 
0 otherwise 


where Z^ = boundary between two consecutive regions 

Zq = 2500 kilometers (upper boundary of the Jacchia- 
Roberts model) 

Z^ = 0 (ground) 

In particular, H^q(Z) = h(Z-Zj^) and = h(Z 2 ~Z). 

Then a model vertical profile can be written as 

4 

F^2i) =2. F^- (Z-) (4-22) 

i'=l 


In this equation, the F^(Z) (i=l-»-4) for temperature, pressure, 

and density are defined in Equations (3-3) through (3-7) . 

For convenience, these equations are repeated here: 

Temperature 


F,(z.) = T„-(Tco-Tf)e 
ECz) = Q] 


-ECZ) 


B Z -Zf 1 ?a t zf 


Tco-V 


(4-23) 




F^tz) = +Zq„(Z-z®,)’' U=z~^) (4- 


24) 


71-1 


•B 


f: = t: 


6 
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Density 


F,(z) =-/«( fjR /f>J 


Pressure 


F,(Z)= - In {f'WpJ 

where 

PjR = R'T ( z > ipok „) 

R T F7R / fCJ ^ 2 ) z 4 loo Km ) 


The F^(Z) (i=2^4) for pressure and density take the same form 

as these for the temperature. The coefficients C. (n=l->-N.; 

B B ^ 

i=2->4) and the boundary values T^, T^, and , can be re- 
adjusted through least-squares fitting procedures. A model 
vertical profile defined by these N unknown parameters is^ 
fitted to a reference profile of observation data defined in 
Section 4. 3. 2.1. These unknown parameters will be collec- 
tively denoted as S , an N-dimensional vector, and F(Z) will 
be written as F(Z,S) to explicitly indicate the dependence 
of F(Z) on S. 
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The reference observation profile is given as 


— 0R .t/? ) I m- M j 


where 5 represents an atmospheric parameter and M is the total 
number of meteorological measurements included in the set. 

The quantities 0^^ and t^^ denote a referencre latitude and a 
reference local solar time, respectively. In the following, 
an abbreviated notation 5(2^^) will often be used in place of 

The loss function to be minimized is defined as 


M 


Q(s)= 21 [F(z„,s) - 4(z„)]' 




(4-27) 


where w_ is the weighting factor of the mth observation. 

m 3 3 

Using matrix notation, this equation can be written as 


Q(s) = [ FCV-i fv [F(s)-i'] 


(4-28) 


where F and Z are M-dimensional column vectors with F(Zj^,S) 
and C(Zj^) as their mth components, respectively, and W is an 
MxM diagonal matrix with w^^ as its mth diagonal element. 

It should be noted that, in Equation (4-28) , ^ should be re- 
placed by ■*^he case of pressure or density. 
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In general, F(S) is not a linear function of S. The linear- 
ized loss function, 0^(5) can be obtained using the Taylor 
expansion of F(S) as follows: 

6 — ( '-MxM matrix 

— _ _ *-3<» (4-29 

AS = S-Se, 




Q^(AS) = (BAS-2i§ ) yf (BAS-A^) (4-30) 


where 

— § - pc So) 


The value of AS that minimizes Q (AS) is given by 

AS = (4-31) 


From this AS, S can be determined as follows: 


S = So + A S 


(4-32) 
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When F(S) is not a linear function of S, this process will 
repeat with S as a new S ^ • If F(S) is a linear function of 
S, i.e., F(S) = Fq + BS, then S which minimizes Q(S) (Equa- 
tion (4-28)) is given by 

S = ( b'^'W B ■if (4-33) 


where 


= € - F. 

The quantity B is calculated as follows: 


B = ( ) 

(4-34) 


Using Equation (4-22 ) , B^^ can be written in terms of F^(Z,S), 

i.e.. 




(4-35) 


4-27 



The solve-for parameter array, S, is defined as follows: 


Cj = (T,®, (??, T„) 

C’l = C Cit , Oz, Ci^i- ) {’t = 2 «v 4 ) 


(4-36) 


Temperature expressions given in Equations (4-23) and (4-24) 
can be written more explicitly as follows: 


F,(Z) = T„-(T„-T,®) 


(4-37) 


B(z) = cj; 


B z-zf Ra-tz! 


B 


Ra't Z 


(4-38) 




Fj(Z) = F,° + Z Cz„ (z-zf)” 


n=i 


(4-39) 


N, 


F,(z) = F,® 




M, 


B 


+ Z, Cj^ ( z - zS ) 

TI=I 


(4-40) 
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IVi 

F^(z) = Ff +Z CjnCzl-z®)" 

n=i 

Z, C3«(z® -z®)" 

n=i 

Q “M 

C4„(z-z!) 


(4-41) 


where 




From these temperature expressions, the partial derivatives 
of temperature with respect to S can easily be obtained, i.e.. 


*9 Ft 

IEl 

'dSz 

2h 


(f + E(z))e 


-ECZ) 


Rx't'L 


i -( 1-h E(z)) e 


-ECZ) 


?Sr, 


0 ( 'n > 3 ) 


(4-42) 
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(4-43) 


lEi - r 


( n =(„3) 


25 

'^^zn 


(Z-Zf)"' (■n=(~Mj) 


9Fx 


( t = 3,4-) 


25 



( ■n = t~3) 


15 

dFs 

^ (^3n 

^r, 


Czf-zf)" 

(z-zf/ 


1-Mx) 




(4-44) 


'^c,n 



C 71=1^3 ) 


2£l 



(t)= t-Ni) 


(4-45) 
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11 

(zt-zlf 

(11=1^ ) 


'SR, 

OQ 

1 

C'n- 

(4-45) 
(Cont ' d) 


In these equations, denotes the Kronecker delta. 

In the coinputation of B , Z will assiame the value of Z . It 
^ mn m 

should be noted that the (3F^/8S^) (i=2->4) are independent of 

S. The partial derivative OF^/3Cj^) is the only partial der- 
ivative that is dependent on S. When the lower atmosphere is 
readjusted independently of the upper atmosphere, the problem 
becomes a linear one and the solution will be given by Equa- 
tion (4-33) . 

For densities and pressures, a similar set of partial deriv- 
atives can be derived. In this case, is redefined as 

c7 = ( T,®, T», ) (4-46) 


where Pq is the density at Z = 90 kilometers. The variation 
of pQ will give an overall rescaling of densities and pres- 
sures. The Fj^(Z) for the density and pressure is given in 
Equations (4-25) and (4-26) , respectively, i.e. , 


Density; Ff (Z) = - ) 


Pressure : 


F. Cz) = - 
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Density : 


(4-47) 


_ __ 1 / ^ fjR 9 fjR '9 fjR '^fjR ^ 

'ic, Pjr\ 9T,“ ’ ’ ST*,’ -afi, y 

Pressure : 

^ a(j® ^ 9To> ^ 9/^./ 


(4-48) 


Since p^ is an overall scaling factor, (9 Pjj^/9Pq) and 
OPjr/SPq) are equal to Pjr/Pq and Pjr/Pq^ respectively. 

The remaining partial derivatives are the same as those given 
in Equations (4-42) through (4-45) , with the following excep- 
tions : 


'dC, ( ' Z = zf 


(for i=2-^4) 


(4-49) 


The right-hand side of this equation is given by Equations 
(4-47) and (4-48) for density and pressure, respectively. If 
the lower atmosphere is readjusted independently of the upper 

g 

atmosphere, then will become an adjustable parameter and 
the same form of partial derivatives derived for the temper- 
ature of the lower atmosphere can be used. 

The computation of (9F^/9C2^) defined in Equations (4-4 7) and 
(4-48) can be carried out analytically, but it involves 
lengthy algebraic expressions. It may be advantageous to 
calculate these partial derivatives numerically. 

If only the lower atmosphere is readjusted independently of 
the upper atmosphere, the upper and lower parts of the atmos- 
phere can be linked using an interpolation scheme. An 
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alternative method would be partially readjusting the Jacchia- 

Roberts model using a new set of parameters associated with 

the lower boundary of the model, i.e., Zq, Pq, and Tq (the 

position, density, and temperature of the last temperature 

minimum near 90 kilometers) . These parameters can easily be 

found from the fitted model lower atmosphere. The atmosphere 

above Z = Z will then be described by the Jacchia-Roberts 
o 

model, and the newly determined lower atmospheric model will 
describe the region below Z = Z^. 

4.3.3 COMPUTATION OF ATMOSPHERIC PARAMETERS 

The final readjusted model vertical profile of an atmospheric 
parameter is given by Equation (4-22) , i.e., 

4 

F(z) = X 


This profile corresponds to a reference latitude, 6^^, and a 
reference local solar time, t^^. Thus, F(Z) can be written 
as F(Z,0j^,tj^) . 

It is assumed that the position of the Shuttle is given by 
(Z ,6 ,(J) ,t). First, the longitude, ()) , and the time, t, 

& o o o 

must be converted into a local solar time, t . Then, 

s 

F(Z ,6 ,t ) can be computed from the following equations, ' 


using F(Z^,0„,t_) and the diurnal, semidiurnal, and latitud- 

S K i\ 

inal variations derived in Section 4. 3. 2.1: 

F(z,,e5,t) = F(Zs,es,t,) oCt^.t,) 


FCZs.ej.t,) = F(z,,e,,t,)+5(z,,0R,t,)(e,-e») (4-5 d 


where D and g are described in Section 4. 3. 2.1. 
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4 . 4 OUTPUT 


4.4.1 ATMOSPHERIC PARAMETER USER FILE 

The Atmospheric Parameter User File will consist of the 
input Shuttle trajectory (altitude, latitude, longitude, 
and time) together with the calculated values of the 
atmospheric parameters (temperature, pressure, density, and 
wind components) as determined by the models for each point 
on the trajectory. The file header record will indicate 
whether the default or fitted models were used in the calcu- 
lations. Table 4-5 contains the format of this file. 

4.4.2 ATMOSPHERIC PARAMETER OUTPUT REPORT 

The Atmospheric Parameter Output Report will be a printed 
record of the information contained in the user file. It 
will be produced at the same time that the user file is pro- 
duced. The contents of this report are the same as those 
displayed in Table 4-5. 

4.4.3 LOCAL ATMOSPHERIC PROFILE REPORT 

A report consisting of one or more vertical profiles of temp- 
erature, pressure, density, and wind components will be an 
output option available to the user. These profiles will be 
determined from the permanent or working files and will be 
reported at the latitude, longitude, and time the user has 
chosen. This choice will be limited in latitude to a re- 
gion centered on 30 degrees north latitude because of the 
design of the m.odels. The local atmospheric profile report 
will serve as a tool in judging the quality of the system's 
output. The information shown in Table 4-5 will be listed 
in the report. 
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Table 4-5. Contents of Atmospheric Parameter User File" 


WORD 

SYMBOL 

DESCRIPTION 

UNITS 

1 

LAT 

Latitude 

Degrees, + N 

2 

LON 

Longitude 

Degrees, + E 
to 360 

3 

ALT 

Geodetic altitude 

Kilometers 

4 

T 

Time (GMT) 

Hours, minutes, 
seconds 

5 

TE 

Temperature 

Degrees K 

6 

PR 

Pressure 

Millibars 

7 

D 

Density 

Grams/cm 

8 

U 

North-south wind 
velocity 

Meters/second , 
+ N 

9 

V 

East-west wind 
velocity 

Meters/second , 
+ E 

10 

WS 

Total horizontal 
wind speed 

Meters/second 

11 

WD 

Direction wind is 
coming from rela- 
tive to true 
north (0 degrees) , 
increasing posi- 
tively clockwise 

Degrees 

12 

TEU 

Estimated uncertain- 
ty in temperature 

Degrees K 

13 

PRU 

Estimated uncertain- 
ty in pressure 

Millibars 

14 

DU 

Estimated uncertain- 
ty in density 

Grams/cm^ 

15 

UU 

Estimated uncertain- 
ty in north-south 
wind 

Meters/second 

16 

VU 

Estimated uncertain- 
ty in east-west 
wind 

Meters/second 


All quantities are real variables one word in length. A 
single record is 16 words long. One record will be pro- 
duced for each point on the input ENTREE trajectory. 
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4.4.4 MODEL ADJUSTMENT REPORT 


After the working files have been built, a model adjustment 
report will be printed. This report will contain statistical 
inforroation concerning the fitting of the models, such as 
residuals, dispersion of the observations in time and space, 
and the number of observations on which the model adjustment 
was based. 

4.5 ANALYSIS OF SOFTWARE REQUIREMENTS 

This section consists of a set of data flow diagrams 
(Figures 4-1 through 4-9) and process specifications which 
form a structured analysis of LAIRS. The LAIRS software 
will be designed in accordance with this analysis . In addi- 
tion, the software development will be undertaken in two 
distinct phases. In Phase 1, the model parameter calcula- 
tion software and a set of default (permanent) files will 
be developed and delivered, allowing the user to test the 
system in its default mode. In Phase 2, the model adjust- 
ment software will be developed. Following the delivery 
of this software, the user will be able to utilize all 
aspects of LAIRS, including the processing of meteorological 
input data. 
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Figure 4-1. Context Diagram for LAIRS 











Figure 4-2. LAIRS Overall Design 
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Figure 4-3. Process 1 Details 
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Figure 4-5. Process 2.3 Details 
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Figure 4-6. Process 3 Details 



Figure 4-7. Process 3.4 Details 




Figure 4-8. Process 3 
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4.3 Details 
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Figure 4-9. Process 3.4.5 Details 
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DATA DICTIONARY 


COEFFICIENTS A set of coefficients required by the atmospheric 

model for determining the altitude profiles of T^ P, 
and as well as coefficients for determining the 
diurnal and semidiurnal (D&SD) variations in these 
quantities and in U and V. Three varieties of 
coefficients exist in LAIRS ; 

DEFAULT COEFFICIENTS - Coefficients determined from 
the default models and stored on permanent files. 

ADJUSTED COEFFICIENTS - Coefficients obtained by 
adjusting the default coefficients on the basis of 
the input meteorological data. These coefficients 
are stored on working files. When no meteorological 
data is available, the set of default coefficients 
corresponding to the date of the Shuttle's reentry 
are stored on the working files after Jacchia- 
Roberts and CIRA merging has been performed. 

NEW COEFFICIENTS - Coefficients resulting from perma- 
nent file maintenance. Special input is required to 
produce these coefficients, which are placed on new 
permanent files. It is not expected that these coeffi- 
cients will be required on each run. 

DATA The special input required by the permanent file 

maintenance option. 

DATE The date of the Shuttle reentry input to LAIRS. 

HEADER INFORMATION Information on reentry profile, models used, etc., 

written to the header of the atmospheric parameter 
output report and the atmospheric parameter user file. 
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DATA DICTIONARY 



F F 

^ 10 . 7 ' 10.7 


Jacchi a- Roberts model constants (see symbol definitions) 


LATITUDE GRADIENTS 


METEOROLOGICAL PROFILES 
MODEL CONSTANTS 


The gradients of T, P, S, U, and V with respect to 
latitude. These are tabulated for various altitudes 
and stored on the permanent and working files. 

5q(Z, (J), t^) - See PARAMETER SET 

Jaechia-Roberts constants 


OPTIONS A set of input options to be supplied by the user on 

keyword cards, including information such as: 

Which model will be used for parameter calculation 
(default, adjusted, or interpolation) 

Which type of input trajectory is to be processed 
(ENTREE Shuttle trajectory, local atmospheric 
profile "trajectory", or both) 

Which wind model will be used (default, adjusted, 
interpolation, or analytic) 

Whether permanent file maintenance is desired 
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DATA DICTIONARY 


PARAMETER SET A set of atmospheric parameters T, P, P, WS , and WD 

(or T, P, P, U, and V) denoted by the symbol These 

are either observed or calculated: 

5q(Z, 0^, (j)^, T - The ith set of observed param- 

eters from the meteorological observations file. 

It is understood that Z covers the entire altitude 
range of this set of profiles. 

5 (Z , 0 , , X ) - The parameters calculated 

C S S 5 S 

from the atmospheric model for a point on the 
Shuttle's trajectory (subscript s) . Here Z^ 
represents only one point, not a profile. 

Note: The subscripts on the arguments of ^ are very 

important. Besides i (ith profile) and s 
(Shuttle trajectory point) , it is possible to 
have R as a subscript. This denotes the refer- 
ence point of the model, 0„ and tj^ ((}> and x are 
transformed to t) , or, in the case of Z, the 
set of reference altitudes used in computing 
the diurnal and semidiurnal coefficients and 
latitude gradients. Thus ^ (Z , 0 , t„) is a 

set of profiles at the reference latitude and 
solar time and computed at the reference alti- 
tudes . 
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DATA DICTIONARY 


SHUTTLE TRAJECTORY The set of points Z„, 0 , (p„, x„ comprising the 

». , I . .. I. . . . S 5 S S 

Shuttle trajectory residing on the ENTREE output 
trajectory . 

STATISTICS Residuals, standard deviations, observation disper- 

sion information, and other information relating 
to the model adjustment process 

WIND U and V Tabulations of the wind components U and V for various 

altitudes. These are stored on the permanent and work- 
ing files as functions of Z, 0^^, and t^^. 



LIST OF SYMBOLS 


Symbol 


Description 


Z 

6 

<P 

T 

t 

T 

P 

P 

WS 

WD 

U 

V 

c 



10 . 

10.7 


altitude 
latitude 
longitude 
time (.GMT) 

solar hour angle (local solar 
time) 

temperature 

pressure 

density 

wind speed 

wind direction 

eastward wind component 

northward wind component 

parameter set consisting of T, 

P, p, WS, and WD, or of T, P, 
p, U, and V 

- observed 

- computed 

geomagnetic index 

10.7-centimeter solar flux 

81-day average of 10.7-centimeter 
flux 
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PROCESS SPECIFICATION 


PROCESS 1; EXECUTIVE CONTROL 

1.1 READ the input keyword card deck. 

STORE the option flags read from cards 

in COMMON. 

DETERMINE from the option flags whether the 

user has requested permanent file 
maintenance. 

IF permanent file maintenance has been 

requested^ then 

PASS to Process 1.2. On return, 
proceed. 

ELSE no permanent file maintenance has 

been requested, so 

PASS to Process 2 for the building 
of working files. On return, 
from Process 2, 

PASS to Process 3 for parameter 

calculation. On return from 
Process 3, 

STOP 


1.2 

READ 

the special input supplied by the 
user for permanent file updates. 


UPDATE 

the permanent file. 


WRITE 

the new permanent files. 


RETURN 

to Process 1.1. 
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PROCESS 2: MODEL ADJUSTMENT 


2.1 DETERMINE 
ELSE 

2.2 READ 

WRITE 

ADJUST 

WRITE 

2.3 

2.3.1 DETERMINE 
IF 

ELSE 

DETERMINE 


from the option flags whether the 
default or adjusted model has been 
requested. 

the default model has been requested, 
then 

PASS to Process 2.2, 

the adjusted model has been requested, 
so 

PASS to Process 2.3. 

from the permanent files the T, P, p, 

coefficients, the latitude gra dients , 

the diurnal and semidiurnal coefficients, 
the U and V profiles, and the_Jacchia- 
Roberts constants K^, F^^q ^ 

corresponding to the data passed 
from Process 1. 

the retrieved latitude gradients, 

diurnal and semidiurnal coefficients, 
and wind U and V profiles to the 
working files . 

the Jacchia-Roberts model coefficients 
(which depend on K^, F^q and F^q 

to match the CIRA model coefficients. 

the (adjusted) default model coeffi- 
cients to the working files. 


if diurnal, semidiurnal, and latitude 
variations are to be modeled. 

no such variations are to be modeled, 
PASS to Process 2.3.4. 

such variations are to be included, 
so 

if more than one profile exists. 
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IF 


ELSE 


2.3.2 READ 


IF 


ELSE 


only one profile exists, 

PASS to Process 2.3.2. 

more than one profile exists, so 

READ the altitude (Z) values of each 
profile. 

EDIT the profiles with too few Z 
values . 

CONSTRUCT a set of reference Z values 
(Z„) from those read. 

I\ 

INTERPOLATE each profile to obtain 
-C (Zr, (})^, T^) . 

TRANSFORM the units of the inter- 
polated profiles to the standard 

model units (tfc., T.->t.). 

^ 1 ^ 1 1 

READ the permanent latitude gradient 
file for each profile parameter. 

TRANSLATE ? (Z„, 6,, t. ) to 
C (Zj^, t^) using the retrieved 

default gradients. 

COLLECT statistics on Z„ and 
C (Zj^, t^) and an original pro- 
file Cq (Zj^, (j)^, t^) , and store 

in COMMON. 

PASS to Process 2.3.2. ~~ 

the permanent files to obtain the 
default diurnal and semidiurnal 
(D&SD) coefficients. 

only one set of profiles exists, 

WRITE the default D&SD coefficients 
to the working file. 

PASS to Process 2.3,3. 

more than one set of profiles exist, 
so 

SOLVE FOR the D&SD coefficients by 
the method of Section 3.2. 
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2.3.3 READ 
IF 


ELSE 


PASS 

2.3.4 READ 


IF 


CHECK the new values against the 
default values, and set a warning 
flag if they are too discrepant. 

WRITE the new D&SD coefficients to 
the working file. 

COLLECT and store statistics on D&SD 
fitting. 

CONSTRUCT C(Zj^, 6^^, t^). 

PASS C(Zj^, 6j^, t^) and 9^, t^) 

to Process 2.3.3. 

the permanent latitude "gradient f il¥T.' 

only one set of profiles exists^ 

WRITE the default gradients to the 
working file. 

more than one set of profiles exists, 
so 

CALCULATE the new latitude gradients 
by the method of Section 4.3. 

CHECK the new gradients against the 
old and set a warning flag if they 
are too discrepant. 

WRITE the new gradients to the work- 
ing file. 

to Process 2.3.4 

the original meteorological observa- 
tion files to obtain the sets of 
profiles Cq (Z^, 9^^, (pj^, t^) . 

D&SD and latitude gradients have not 
been requested, 

MERGE these profiles to obtair^the 
single reference profile set 

5 (Z . , 9„, t^) - assume D&SD and latitude 

gradients are zero. 
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ELSE 


IF 


ELSE 

PASS 


2.3.5 READ 


DETERMINE 


IF 


ELSE 

SOLVE FOR 


D&SD and latitude gradients have been 
modeled and are available in COMMON, 
so 

TRANSLATE these profiles to obtain 
^(Z., 0 , t„) using the D&SD coeffi- 

cients and the latitude gradients 
(default values are used when there 
is only one profile set) . 

winds have been requested, 

WRITE the U and V profiles from the 
5(Z., 0 , t_) set to the working 

1 K i\ 

files . 

do not write to wind files. 

T, P, P from the C(Z., 0„, t„) set 

IKK 

to Process 2.4. 

the permanent files to obtain the 
default model coefficients T, P, P 
for all segments of the atmosphere. 

if any default points must be added 
to the reference profile, either 
because of user request or to make 
the profile complete. 

any points must be added, then 

CALCULATE the needed points from 
the default model coefficients and, 
if the Jacchia-Roberts section is 
involved, from the K^, ^10. -j , and 

F^q ^ values read from the permanent 
Jacchia-Roberts files . 

ADD these points to the reference 
profile. 

the profile is complete as is. 

the adjustments to the coefficients 
using the iterative scheme discussed 
in Section 4.3. 
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WRITE 


the adjusted model coefficients to 
the working file. 


COLLECT 

WRITE 

RETURN 


statistics on the model fitting 
process . 

the statistics to the model adjust- 
ment statistics report. 

to Process 2.1. 
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PROCESS 3: ATMOSPHERIC PARAJ4ETER GENERATION 


3 . 1 DETERMINE 


IF 


ELSE 


3.2 READ 
PERFORM 
PASS 

3.3 CALCULATE 
PASS 

3.4 

3.4.1 DETERMINE 


IF 


from the option flags whether the 
parameters to be calculated are for 
an ENTREE trajectory or are for a 
Local Atmospheric Profile (more than 
one LAP may be requested; if so, this 
begins an outer loop) . 

an ENTREE profile is requested, then 

WRITE the Atmospheric Parameter Report 
header. 

PASS to Process 3.2 

a Local Atmospheric Profile is 
desired, so 

WRITE the LAP report header. 

PASS to Process 3.3. 

a trajectory point from the ENTREE 
trajectory file. 

any needed transformation of units 
for this point. 

to Process 3.4 for parameter calcula- 
tion for this point. 

a "trajectory" point for the vertical 
Local Atmospheric Profile. 

to Process 3.4 for parameter calcula- 
tion for this point. 


from the trajectory point passed 
from Process 3.2 (ENTREE) or' 373’ 

(LAP) whether the upper or lower 
atmospheric model is required (check 
value of Z) . 

Z lies within the bounds of the upper 
atmospheric model, then 

PASS to Process 3.4.2 (upper atmos- 
phere control) . 
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ELSE 


3.4.2 DETERMINE 


IF 


ELSE 


PASS 


Z lies within the bounds of the 
lower atmospheric model , so 

PASS to Process 3.4.4 (lower atmos- 
pheric control) . 

if this is the first time the upper 
atmospheric model has been called 
for this trajectory. 

this is the first call, the model 
needs to be initialized, so 

READ the working upper atmospheric 
files . 

STORE the retrieval coefficients in 
COMMON . 

CALCULATE the model constants 
and store. 

the coefficients have already been 
retrieved, so proceed. 

to Process 3.4.3 for calculation. 


3.4.3 DETERMINE if the interpolation model has been 

requested. 

IF the interpolation model has been 

requested, then (3. 4. 3. 3) 

SEARCH for the meteorological profile 
nearest this point, and retrieve 9 
and T for that profile. 

INTERPOLATE this profile to find 
T, P, p, WS, and WD for this altitude. 

TRANSFORM units of parameters to model 
standard units. 

TRANSLATE the retrieved parameter 
set in hour angle and latitude to 
the position of the trajectory point 
using D&SD coefficients and latitude 
gradients (if requested) . 
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ELSE 


ELSE 

3.4.4 DETERMINE 

IL 

ELSE 

3.4.5 DETERMINE 
IF 


the Jacchia-Roberts model will be 
used, so (3. 4. 3.1) 

CALCULATE the values of T, P, p for 
this trajectory point using the J-R 
equations, the retrieved model coeffi 
cients, and constants. 

Winds U and V have been requested, 
then (3. 4. 3. 2) 

RETRIEVE U and V from working wind 
file for this point. 

No wind model has been requested, 
so , 

PROCEED to Process 3.4.6 for 
output. 


if this is the first time this seg- 
ment of the lower atmospheric model 
has been entered on this trajectory. 

this is the first time, the model 
needs to be initialized so, 

READ the working files for this 
model segment. 

STORE the retrieved coefficients and 
gradients in COMMON . 

the coefficients for this segment 
have already been retrieved, so 

PROCEED to 3.4.5 for parameter 
calculation. 

if the interpolation model has been 
requested. 

the interpolation model has been 
requested, then (3. 4. 5. 4) 

SEARCH for the meteorological profile 
nearest this trajectory point. 

RETRIEVE 6 and x for this profile. 

INTERPOLATE this profile to find 
T, P, C, WS, and WD for this altitude 
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TRANSFORM units of parameters to 
model standard units. 


ELSE 


DETERMINE 

IF 

IF 


ELSE 


TRANSLATE the retrieved parameter 
set in hour angle and latitude 
(t and 6) to the position of the 
trajectory point, using D&SD coeffi- 
cients and latitude gradients . 

the polynomial models have been 
requested, so 

CALCULATE the values of T, P, P for 
the altitude from the retrieved model 
coefficients . 

TRANSLATE the calculated parameters 
T, P, p to the t and 9 of the traj- 
ectory point using the retrieved D&SD 
coefficients and latitude gradients 
(if requested) . 

if wind modeling has been requested. 

no wind requested, 

PROCEED to 3.4.6 for output. 

winds from working files requested, 
then, 

RETRIEVE U, V, and latitude gradients 
for this point from working file 
(3. 4. 5. 2) . 

!I;raNSLATE U and V to the 0 and x 
of this trajectory point, using 
retrieved D&SD coefficients and 
latitude gradients. 

analytic wind model has been 
requested, so 

CALCULATE U and V from the wind 
equation, using T, P, and p pre- 
viously calculated and gradients 
previously retrieved from working 
file. 

PROCEED to 3.4.6 for output. 
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3.4.6 WRITE 


DETERMINE 


IF 


ELSE 


the values of the trajectory point 
(Z, 6, (j> , t) and the corresponding 
calculated parameters (T, P, p, and 
horizontal wind compon^i^ts ) to the 
Atmospheric Parameter Report^_or the 
Local Atmospheric Profile !Report._ 
Also write to the user file i f this 
is an ENTREE, not a LAP, trajectory^ 


if the point just processed is the 
last point on the LAP or ENTREE 
trajectory. 

this was the last point, this parti- 
cular profile is completed, so 

RETURN control to Process 3.1 
t atmospheric model executive) to 
begin processing the next requested 
ENTREE or LAP trajectory. 

more points on this trajectory 
remain to be processed, so 

RETURN control to Process 3.2 (if 
this is an ENTREE trajectory) or 
Process 3.3 (if this is a LAP) to 
process the next trajectory point. 
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4 . 6 SOFTWARE CONSTRAINTS 


4.6.1 EXECUTION ENVIRONMENT 

The Langley Atmospheric Information Retrieval System will 
be designed to operate on the CDC 6000 and Cyber series 
computers at the Langley Research Center at Hampton, Virgin! 
The upper limit on the main memory core requirements should 
be approximately 200K. 

The language of LAIRS will be high-level FORTRAN consistent 
with LaRC coding standards. No assembler code will be 
written. No routines will contain multiple entry or exit 
points, nor will they contain assigned variables. The use 
of GO TO statements will be held to a minimum. All COMMON 
blocks will be as small and as specialized as possible; 
there will be no use of blank, global COMMON. Subroutines 
of unreasonable length will not be permitted. 

4.6.2 MAINTENANCE CONSIDERATIONS 

The LAIRS development will utilize modular software through- 
out, so that new capabilities can be easily added to the sys 
tern. An attempt will be made to foresee possible areas of 
future enhancement and to provide appropriate dummy subrou- 
tines or, if dummy routines are not feasible, easily alter- 
able code. 

The model adjustment software will be designed to allow the 
user to rebuild the permanent files should superior default 
data become available. The Jacchia-Roberts , K , F,n , and 
Fj^q ^ files will require periodic updating, which may be 
accomplished by software available from CSC. 

Consistent subroutine and COMMON block documentation will 
be stressed at all stages of development. This documenta- 
tion, combined with the modular structure of the system 
and the coding standards elaborated upon in Section 4.6.1, 
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will insure that LAIRS is an easily maintainable, upgradable 
system. 
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